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Atmospheric  Attenuation 
of  Millimeter  and  Submillimeter  Waves: 

Models  and  Computer  Code 


1  INTROIUKTIOiN 

The  increasing  use  of  mici'cnvavos  in  the  past,  few  years  by  DOD  for  determin¬ 
ation  of  geophysical  parameters,  high  data  rate  communication  channels,  weapons 
guidance  and  target  acquisition  lias  emphasized  the  need  for  better  modelling  and 
more  accurate  calculations  of  attenunviyn  through  clear  and  iiyd rometeor  laden 
atmospheres.  To  fulfill  this  need,  the  use  of  millimeter  and  submillimeter  wave¬ 
lengths  show  many  advantages  over  the  use  of  microwaves: 

a.  better  sensitivity  to  geophysical  parameters 

b.  larger  communication  band  widths 

c.  more  secure  communications 

d.  narrower  beam  widths 

o,  good  angular  resolution 

f,  imagery  capability 

One  disadvantage  is  atmospheric  attenuation  caused  by  gaseous  absorption 
and  extinction  by  hydrometeors. 

The  objective  of  this  report  is  to  model  and  calculate  atmospheric  transmis¬ 
sions  and  atmospheric  attenuation  of  radio  frequency  waves  in  the  spectral  region 
of  l -1000  GHz,  that  is,  30-.  03  cm  wavelengths  or  .30-33  cm  *  (wave  numbers). 

(Received  for  publication  15  October  1979) 
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}  liis  ropori  lollows  in  i li<*  spiril  « >i*  l li<*  nn^iiub  atnu>splioric  transmission  niodols 

by  MrC’latrlioy  rt  alf  1  wlnoh  woro  d«'v«*mptcl  into  flu*  ),(A\TI\AN  transmission 
2  * 
oodof  an  oj>L-rai  ionally  dirootrd  oompntt'r  uodt*  1\>i*  user  iniplomonl  at  ion. 

Tlio  modt'ls  usod  in  tins  study  arc  tho  Ch  ar  Aimosphiwt*  Modrl  ^ivo.n  in 
Valley,  tin-  log  ami  Cioud  Models  from  Silverman  ami  the  research  of  the 
authors  (lexis  In  Mason’’  and  Borovikov*’  v.ito  extensively  utilized),  and  llu-  Urdu 
Models  which  were  generated  li-ym  Hie  research  of  the  authors.  The  rainfall  rate 
considered  are  from  0.2a  nun, 'hr  to  laO  mnt/hr.  The  drops  i/o  distribution  of 
these  rainlall  rates  is  determined  from  the  Marshall-Palmer  distribution. 

The  use  ol  these  models  should  not  be  Interpreted  as  if  these  are  the  only 
models  ot  use  in  propagation  studies.  The  models  were  chosen  because  thev  are 
representai  ive  ol  niirilatitude  conditions.  The  computer  program  given  in 
Appendix  P  will  allow  eaeh.  researcher  to  calculate  transmittance  and/or  attenua¬ 
tion  lor  any  horizontally  stratified  model  which  represents  real  world  conditions. 

Table  1  lists  the  delinitions  for  the  spectral  region  c  the  electromagnetic 
spectrum  that  will  be  employed  in  this  report. 


Table  1.  Spectra  Designations 


Region 

Wavelength 

Frequency 

Microwave 

1  m  -1  cm 

.  30G1IZ-30C11?. 

Millimeter 

1  cm  -0.  1  cm 

30-300 

Submillimeter 

0.  1  cni-0.  01  cm 

300-3000 

The  codes  described  in  (his  report  are  available  from  National  Climatic  Center 
Federal  Building,  Asheville,  NC  28801  for  a  service  charge. 

1.  MeClatchey,  H.A.,  Venn,  11. W.,  Selby,  J.,  Volz,  F.E.,  and  Caring,  J.S. 

(1972)  Optical  Properties  of  the  Atmosphere  (3rd  lid.)  A  FCR1  ,-72-0497 
A  D 

2.  Selby,  J.R.A.  and  MeClatchey,  H.A.  (1973)  Atmospheric  Transmission  From 

0.  2a  to  j 3.  3  um ;  Computer  C  ode  ROWTRAN  2.  A !•'(' K I 7  2 -07 43 
AD - 

3.  Va^Uty.^S.  B.  (Fditor)  (19ti5)  Handbook  of  Geopliysics  and  Space  Knvi ronments, 

4.  Silverman,  B.  A ,  and  Sprague,  E,  D.  (1970)  Airborne  Measurements  of 

In-Cloud  Visibility,  National  Conference  on  Weather  Modification  of  the 
American  Meteorological  Society,  April  ti-9,  Santa  Barbara,  California. 

o.  Mason,  B,  J.  (1971)  I  he  Physios  of  Clouds,  Clarendon  Press,  Oxford. 

G.  Borovikov,  A.M.,  Kingian,  A.K.H.  and  others  ( 1063)  Cloud  physics,  P.S. 
Department  of  Commerce,  Office  of  Technical  Services. 
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Tlic  term  radio  frequency  encompasses  the  min  nwavc,  milluucl.  r  and  suhunlli- 
molt'i  frequency  bands. 

In  the  1-1000  GH /  frequency  range  considered,  the  clear  atmosphere  attenua¬ 
tion  is  tiue  primarily  to  the  rotational  absurptkm  lines  of  water  vapor  (ll.,ti)  and 
oxygen  (O.,).  It  one  were  mainly  interested  in  communication  channels,  only  water 
vapor  and  oxygen  are  needed  to  account  for  clear  air  attenuation.  However,  the 
applieat  ion  of  this  research  in  addition  to  communication  is  the  identification  of 
geophysical  parameters,  ground  target  acquisition  and  weapons  guidance.  These 
.applications  require'  the  additional  absorption  lines  of  ozone  (O.,)  carbotl  monoxide 
(CO)  and  nitrous  oxide  (\,,o)  lor  completeness  and  accuracy. 

Since  the  A1G1,  1IITKAN  data  base  and  computational  technique  id'  McClutohey 
el  a)  1  was  originally  genet  ated  by  specialists  in  the  infrared  region,  there  exist 
differences  ill  notation  between  it  and  the  notation  employed  by  the  microwave  spe¬ 
cialists.  In  the  infrared  region,  the  specialists  use  .  ratisinittanee  to  describe 
propagation,  whereas  in  the  microwave  region  the  specialists  use  attenuation  to 
describe  propagation  (see  Appendix  A  for  definitions), 

tine  major  point  of  ibis  research  is  that  systems  designers  lia  -e  to  adjust  their 
ideas  of  completely  deterministic  propagation  models  and  rethink  system  propaga¬ 
tion  as  being  statistical  in  nature,  hi  the  microwave  region  the  atmosphere  is 
essentially  static  (that  is,  at  most,  rain  conditions  have  to  be  considered  for  fre¬ 
quencies  above  10  till:').  As  systems  go  higher  in  frequency  through  the  millimeter 
and  submillimeter  bands,  the  clear  atmosphere  attenuates  due  to  water  vapor, 
oxygen,  ozone  and  other  minor  constituent  gases.  In  addition  to  the  clear  atmos¬ 
phere,  fog,  clouds  and  rain  attenuats  by  belli  scattering  (deflection  of  energy)  and 
absorption  (dielectric  heating  of  water  drops).  It  is  the  statistics  of  the  meteoro¬ 
logical  conditions  of  gaseous  atmospheric  concent  rat  ions,  and  hydrometeor  fog, 
clouds  and  rain  that  have  to  be  considered  when  designing  millimeter  and  sub- 
nbllinieter  systems. 

2.  C.I.KAR  ATMOSPHERE 

The  atmospheric  molecular  absorption  spectra  is  calculated  by  a  computer 
efficient  algoritlun  of  Al’GJ.'s  1I1THAN  Code.  This  new  code  was  advanced  by 
Clough  et  al.  It  is  called  l'ASCOl'-l,  an  acronym  for  l'ast  Atmospheric 

7.  A'lcClatche.V,  11.  A.,  Honed  iet,  VV.  rf. ,  Clough,  S.A.,  Ilureh,  D.  1£.  ,  Cal  fee, 

H.F. ,  Fox,  K. ,  liothman,  I..S.  and  Caring,  J.S.  (1973)  AFOUL  Atmos  - 
pherie  Absorption  Fine  Parameters  Compilation,  AFCI1F-T11  -73 -0090 
7TB - — - - 

d.  Clough,  S.A.  ei  al.  (19H0)  FASCOU-1,  Al’GI.  Technical  Report,  lobe 
published. 

1  1 


.'■dp, nature  Code,  The  111  TUAN  code  is  a  compilation  of  spec!  ml  lino  p,u  'JimU-rj? 

(<1  tin-  signi  firaitl  atmospheric  absorption  inn  s  for  frequencies  from  1  lit'  microuav*' 
through  tho  inf i  arr'il  regions  of  1 1 1  c '  spectrum.  As  previously  stall’d,  the  molecules 
of  interest  in  tin'  1-1000  (Jll/  reg’oll  are  oxygon  (O.,),  water  vapor  (ll.pl),  o/oiv 
(1)  ),  raibon  monoxide  (C'O),  and  nitrous  oxide  The  paper  by  i  lough  al' 

should  lie  consulted  lor  further  i’.iforniat  ion  about  this  program  and  the  theory  of 
clear  aii1  transmission. 


.1.  UVDK0M1  TKOKS 


Haze,  fog,  and  clouds  aic  intimately  related.  A.  is  well-known  the  atmos¬ 
phere  is  never  free  of  particles  (aerosols).  Haze  is  an  aerosol;  it  is  composed 
of  hygroscopic  and  non-hygrosoopio  particles.  The  hygroscopic  particles  act  as 
nuclei  for  condensation  of  water  vapor.  In  essence,  fog,  clouds  an,  .'ms,  precipi¬ 
tation  are  the  meteorological  results  of  nuclei  that  began  as  ha/e.  The  transmit 
tance  and  attenuation  tine  to  hydrometeors  (log,  clouds  and  rail')  is  calculated  by 
the  full  Mie  theory.  The  only  distinction  between  fog  and  clouds  ts  their  distance 
from  the  earth.  Similarly,  the  distinction  between  clouds  (or  fog)  and  rain  is 
almost  non-existent  when  they  are  mixed,  Fui'  conciseness,  the  terms  fog,  cloud 
and  lain  have  to  be  defined.  la  this  report,  the  radar  meteorologists'  definition 
will  be  adhered  to.  Tog  and  clouds  refer  to  a  collection  of  water  drops  with  diam¬ 
eter  less  than  100  gm,  [N.ll.  this  report  uses  radius  instead  of  diameter;  thus  iho 
collection  of  water  drops  with  radii  less  iiinn  at)  pm  or  0.  (la  mm],  Haiti  refers  to 
a  collodion  of  water  drops  with  diameters  up  to  i>.  5  mm  (or  radii  up  to  ,'i .  3  mm). 

It  is  assumed  that  the  hydrometeors  are  homogeneous,  spherical  masses  ,,f 
pure  water,  with  a  density  slightly  less  than  1  g  cm’’,  A  quick  observation  of 
these  definitions  is  that  the  fog  and  cloud  droplets  are  much  smaller  than  the  wave¬ 
lengths  considered  ill  this  report.  11  will  he  shown  later  that  in  this  ease,  the 
llaylcigh  approximation  holds  (x  =  2*r/A  w  1),  and  the  attenuation  of  the  fog  ot 
cloud  is  directly  proportional  to  the  liquid  water  content  (see  Kq.  (21'!.  On  the 

other  hand,  rain  attenuation  requires  the  full  Mie  theory  calculations.  Il  >s  as-  . 

Slimed  that  the  raindrop  spectra  lias  the  distribution  first  advanced  bv  Marshal! 
and  Palmer'  (M-P).  The  choice  of  M-P  distribution  is  mnd.fj&'or  the  following 
reasons:  ”  t 

a.  The  rain  parameter  measured  by  meteorologists  is  the-rainfall  rate  It, 
but  tlie  parameter  of  iriportanoe  in  scattering  calculations*!!?'  dropsizc  dish  ibution. 

V.  .a 


9.  Marshall,  J.S.  and  Palmer.  W.M.K.  (1948)  Tnc  distribution  ol  raindrops  with 


The  M-l>  dist  ribut  ion  n(r)  is  explicitly  related  to  rainfall  rate  H  (see  l'lqs.  (33) 
and  (33)). 

li.  Tlio  M-r  distribution  has  a  negative  exponential  form,  which  is  the  gen  ¬ 
era]  form  or  othei  "accepted"  distributions.  Specifically,  the  same  form  as  the 
]  )ci  rtiiend  p  ar.^  modified  gamma  ftlnelion  distribution  used  in  modelling  cloud  and 
rain  droplet  distributions  where  y  1  (see  Fq.  (lfi)), 

e.  The  amplitude  factor  in  the  M-P  distribution  is  a  const  ift  for  rainfall  lor 
continuous  rain  conditions  (see  Table  li). 

d.  Many  ri  searchers  are  currently  employing  the  M-T  distribution,  so  it  is 
already  accepted  by  the  reseat  eh  community. 

t.Ai.f.n.vnoNs 

4.1  (iii'ai*  A(iiH*glu'io  Iransitiissiim/AKviiiialitm 

In  Ihe  frequency  range  1-1000  till/  the  clear  atmosphere  progresses  from 
100  percent  transparent  to  opaque  as  the  frequency  is  increased.  This  transmis¬ 
sion  and  attenuation  is  present  100  percent  of  the  time  (see  Figures  1,  2,  3  and  1). 
Clough  et  ai°  discusses  the  clear  air  transmittance. 

4.3  Ilydrninrli-tiiV  Attenuation 

Tim  t  i  nnsmim.nce/aUetuiation  of  millimeter  and  submiliimeler  waves  by 
liytit  oineteoru  is  appreciable  for  a  significant  percentage  of  a  year.  Consequently, 
estimates  of  fog,  cloud  and  ram  attenuation  arc  required  for  the  design  of  various 
systems,  Physically,  extinction  (It!  notation  which  is  equivalent  to  attenuation  m 
microwave  notation),  the  removal  of  energy  from  an  electromagnetic  wave,  is  the 
sum  of  two  processes.  The  first  is  absorption  (Q^)  which  arises  from  the  dielec¬ 
tric  heating  of  condensed  water  by  the  incident  electromagnetic  wave,  and  the  sec¬ 
ond  is  scattering  (Q  )  which  arises  from  scattering  of  the  incident  electromagnetic 
wave  into  other  directions  (see  Iiq.  (1)>.  It  should  be  noted  that  the  scattering 
term  includes  radiation  that  is  scattered  into  the  foreward  direction  as  well  as 
into  other  directions.  The  attenuation  due  to  frozen  (ice)  hydrometeors  is  not 
considered  in  this  report,  since  ihe  imaginary  component  of  the  index  of  refraction 
is  small  and  thus  the  attenuation  is  small.  Stated  differently,  natural  ice  ciouds, 
hail  and  snow  are  highly  transparent  to  passive  radiowave  measurements.  It  is 
sale  to  say  that  the  theory  of  calculation  of  scattering  by  hydrometeors  is  well 
known  (see  Figure  5).  However,  the  meteorological  statistics  of  hydrometeors 
are  not  well  known.  This  latter  fact  leads  to  difficulties  in  modelling  as  will  be 
shown  below. 

10,  Deirniendjian,  1J.  (10(19)  Klcctrpmagnctic  Scattering  on  Polydispersions. 

Klsevier  Publishing  Co. 
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Figure  4.  The  Attenuation  Through  a  X  km  Thick  Atmosphere  at  Sea  Level 
(20°C,  1013  mb)  and  4  km  (0°C,  628  mb).  This  figure  relates  attenuation  de¬ 
crease  with  altitude 
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Figure  5.  Attenuation  vs  Frequency  for  Rainfall  Rates,  0.  25,  1.  25,  2.  50, 

5,  00,  12.  50,  25.  50,  100  and  150  mm/hr  for  20°C  temperature  drops  employ 
the  Marshall-Palmer  distribution.  Experimental  data  verifying  theoretical 
calculations  are  noted.  The  experiments  measured  both  attenuation  (dB/km) 
and  the  rainrate  (mm/hr).  Oguchi's  calculations  are  included  for  comparison 


Due  to  the  number  of  references  to  be  included  as  footnotes  on  this  page,  the  read¬ 
er  is  referred  to  the  list  of  references,  page  53. 
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In  most  scattering  calculations  it  is  assumed,  as  in  this  report,  that  the  shape 
of  hydrometeors  is  approximately  spherical  and  their  temperature  is  homogeneous. 
These  two  assumptions  are  true  for  l'og  and  clouds.  'When  water  particles  (fog  and 
clouds)  are  small,  that  is,  .  01-10  ym  in  radius  the  ratio  of  their  surface  to  mass 
is  large,  of  the  order  of  10*’- 10^  respectively.  This  means  that  small  water  par¬ 
ticles  respond  rapidly  to  changes  in  temperature.  Rain  behaves  diffeiently;  it  is 
a  well  known  fact  that  raindrop  shape  depends  on  the  drop  radius  and  the  conditions 
of  fall.  Raindrops  smaller  than  1  mm  in  radius  tend  to  be  spherical.  Drops 
approximately  1.  5  mm  and  larger,  falling  in  stagnant  air,  tend  to  be  flat  or  even 
concave  at  the  base  and  slightly  rounded  at  the  top,  whereas  drops  (--  1.  5  mm), 

5  j  9 

falling  in  turbulent  air,  tend  to  be  ellipsoidal.'  Jones  studied  the  shape  of  rain¬ 
drops  in  naturally  occurring  rain  and  has  classified  them  as  spherical,  oblate  and 
prolate  spheroids,  and  irregular  ellipsoids.  The  large  drops,  he  concluded, 

2  Q 

oscillated  about  a  preferred  shape.  Oguchi  calculated  scattering  from  spheroids. 
A  general  theory  for  other  shapes  does  not  exist.  Rain  has  a  surface-to-mass 
ratio  of  the  order  of  10  (which  means  it  does  not  respond  rapidly  to  changes  in 
temperature)  so  that,  as  the  drop  falls,  heat  transfer  between  the  drop  and  Us 
environment  is  taking  place.  Thus,  the  temperature  of  a  raindrop  is  not  homo¬ 
geneous  because  of  this  heat  transfer  between  the  drop  and  its  environment. 

Other  implicit  rain  assumptions  are  that  the  drops  are  not  electrically  charged, 
and  that  the  drops  do  not  cluster  together.  This  latter  assumption  ignores  the  fact 
that  raindrops  are  not  uniformly  dispersed  in  the  atmosphere,  even  though  there 

is  evidence  to  show  that  raindrops  tend  to  group.  The  usual  assumption  made  is 

2 1 

that  the  drops  are  far  enough  apart  to  have  negligible  interaction  between  them. 

The  dual  assumption  of  spherical  drops  and  homogeneous  temperatures  of  drops 
has  been  made  to  simplify  calculations. 

The  drops  of  liquid  water  in  fog,  clouds  and  rain  absorb,  scatter  and  emit 

thermal  radiation  at  radiowave  frequencies.  The  absorption  and  scattering  prop- 

22 

erties  of  a  single  dielectric  sphere  was  first  advanced  by  Mie.  Fog,  clouds  and 
rain  are  not  composed  of  a  single  size  drop,  but  are  a  composition  of  a  distribution 
of  many  sizes.  IE  this  distribution  of  dropsize  radii  is  known,  then  the  Mie  theory 


19.  Jones,  D.M.A.  (1959)  The  shape  of  raindrops,  J.  Meteor.  16:504-510, 

20.  Oguchi,  T.  (1960)  (1964)  Attenuation  of  electromagnetic  wave  due  to  rain  with 

distorted  raindrops,  J.  Radio  Res.  Lab.  (Japan)  7:467,  Part  II,  J.  Radio 
Rea.  Lab.  JL_1: 19.  - 

21.  Lu'  e,  G.D.  (1968)  Penetrability  of  haze,  fog,  clouds  and  precipitation  by 

radiant  energy  over  the  spectral  range  0.  1  micron  to  10  centimeters, 
NAVWAG  Study  01.  AD  847658,  MaJ. 

22.  Mie,  G.  (1908)  Beitrage  zur  optik  truber  medien,  speziell  kollodaler 

metallosungen,  Ann  der  Phys.  25:377-445. 


may  be  used  to  calculate  the  exact  absorption  and  scattering  properties.  An 
excellent  review  and  application  paper  was  published  by  Gunn  and  Fast,  applying 
the  Mie  Theory  to  the  scattering  of  microwaves  by  liquid  and  solid  hydrumeteors. 
A  limitation  of  the  Gunn  and  East  work  is  that  they  considered  only  those  fre¬ 
quencies  of  use  to  radar. 

In  the  Mie  Theory,  the  following  terms  are  defined  for  the  scattering  of  a 

plane  wave  by  a.  homogeneous,  dielectric  sphere  (pure  water  sphere  in  this  re- 

2 

search).  The  extinction  cross-section  Q  (  (I,  )  is  the  sum  of  the  absorption 
cross-section  Qa  and  the  scattering  cross -section  Q^,  (where  L.  is  the  unit  of 
length) 

QextsQa+Si  <l> 

The  extinction  cross-section  is  defined  such  that  Q  ^  multiplied  by  the  incident 
power  yields  the  total  power  extracted  by  an  obstacle,  cither  by  absorption, 
scattering,  or  both,  from  the  incident  plane  wave. 

_o 

Suppose  that  wr  is  the  incident  power  density  (Watts  1,  )  upon  the  droplet. 

The  extinction  cross-section  is  the  ratio  of  the  total  power  P(watts)  extracted  to 
the  incident  power  density 

^ext  =  W~  (l2)  <2> 

In  the  Mie  Theory  the  scattering  coefficients  are  expressed  in  terms  of  a  dimen¬ 
sionless  size  parameter  x  =  kr,  where  k  =  2 rr/\  is  the  free  space  propagation  con¬ 
stant  and  r  is  the  radius  of  the  sphere 


A.  is  the  wavelength  of  incident  radiation  surrounding  the  particle.  The  physical 
importance  of  the  particle  radius  and  the  wavelength  is  that  these  two  parameters 
determine  the  distribution  or  phase  over  a  particle.  This  is  seen  from  the  incident 
electric  field  Einc  =  exp  [-i(kr  -  ut)]  .  When  the  radius  of  the  particle  is  much  less 
than  a  wavelength,  the  phase  of  the  incident  wave  is  uniform  over  the  particle 
(Rayleigh  scatter).  When  the  size  of  the  particle  is  comparable  to  a  wavelength, 
the  phase  of  the  incident  wave  is  not  uniform  over  the  particle,  and  there  results 
spacial  and  temporal  phase  differences.  These  differences  exhibit  themselves  as 


23.  Gunn,  K.  L.  S.  and  East,  T.W.  R,  (1954)  The  microwave  properties  of 
precipitating  particles.  Quart.  J.  Roy.  Meteor.  Soc.  80:522-545. 
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inter fereneos  in  the  scattered  wave.  The  interferences  in  turn  depend  upon  the 
wavelength  of  incident  radiation,  the  size  (r),  index  of  refraction  (in)  of  tiie  par¬ 
ticle,  and  the  angular  position  (0)  around  the  sear.terer.  It  is  the  relative  size, 
that  is,  ratio  of  the  particle  size  to  wavelength  (r/A),  which  is  the  physically 
important  parameter.  At  lower  relative  size,  that  is,  x  «  1,  Mie  Theory  reduces 
to  the  Rayleigh  approximation  (Figures  C  to  fi).  Specifically,  the  Mie  Theory  is 

a  mathematical  series,  the  first  terms  of  which  arc  equivalent  to  the  Ravleirh 
24  .  o 

expression.  Tims  Mie  Theory  describes  accurately  the  extinction  by  particles 

of  any  relative  size.  The  Mie  extinction  cross -seet ion  Q  ,  and  scatter  cross- 

ext 

section  are  defined  as 


-?T  £  1  <an  r  V 


CD 


%  (2n  +  1)  C|anl2  H  |bn|2) 

n  1 

and  from  1.1  q.  (1)  the  absorption  cross-section  Q  is 

o 

a  oxl  s 

where  a  h  air  the  Alio  coefficients 
Re  denotes  "real  part  of" 

The  Mie  coefficients  arc  given  by 

-in(mx)[x.in(x)]  1  -  ,ii,(x)|mxli  (nix)] 1 

it  - - - - -  ■  1  ‘  - - 

n  jn(mx)[xhj^'<x)] 1  -  h^(x)[mx.iiOnx)l 1 

-j  (x)[mxj  (mx)|  '  -  m2i  (x)[x)  (x)J  ' 

_  1 1 _ n _ n _ n _ 

11  h^(x)[mxjn<m)J  1  -  in2jn(mx)[xh^2^(x)l  1 


(4) 


<  1 1 ) 


<f>) 


(fi) 


Penndorl,  R.  (lfhi:i)  Research  on  Aerosol  Seatierim 
Final  Report,  AFCR I,-H3 -I1HH,  AD 


in  the  Infrared, 
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Kignre  8.  The  Relationship  of  the  Rayleigh  Attenuation  to  the  Temperatures  0, 
20,  and  40°C 


where  in  is  the  complex  index  of  refraction  for  the  particle  (drop) 

'  -  derivative 
2ttr  . 

x  =  — —  size  parameter 

A 

r  -  radius  of  particle 

j  is  n^1  order  spherical  Pcssel  function 
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is  nlh  order  snhuri«-nl  IJanlccl  function  of  the  second  kind, 
n  -  -i 

The  spherical  Bessel  and  IJankel  functions  arc  related  to  the  ordinary  cyclindrieal 
(2) 

functions  j  (x)  and  H  <x)  by 


=  (£r)  J.,+i/2(x) 


hn,(x)  =  (£)1/2*n\/2M 


(7) 


(«) 


Headers  should  be  aware  of  the  differences  in  notation  between  this  report  and 
works  of  Van  de  Hulst,  arid  Deirinendjian.  Van  de  Hulst  employs  an  efficiency 
factor  for  extinction  Qp  ,  defined  as  the  ratio  of  the  extinction  cross-section  to  the 
geometric  area  presented  by  the  drops,  (hat  is 


Q 


ext 


ext 


irr' 


In  Van  de  Hulst's  notation  Eq.  (3)  becomes 


j-2 

2 


23  (2n  +  1)  Re  (an  +  bR) 
n  =  1 


(3-) 


Physically,  the  extinction,  cross-section  efficiency  factor  is  the  fractional  cross- 
sectional  area  of  a  particle  that  acts  on  the  incident  wave.  Similar  definitions  for 
scattering  and  absorption  efficiencies  are  defined  by  Van  de  Hulst.  Deirmendjian 
on  the  other  hand,  uses  K  to  represent  Van  de  Hulst's  efficiency  factors.  In  this 
report  K  -  (m‘  -  l)/(nA  +  2)  (see  Eqs.  (10)  and  (ID). 

In  the  limiting  case  when  r  «  1,  that  is,  x  «  1,  the  power  series  representa¬ 
tion  of  a  and  b  in  terms  of  a.,  b.  and  b9  are  significant  if  higher  terms  than 

n  n  25  1  1  24  e‘  11 

are  neglected,  see  Kerr,  Pennuorf,  or  Van  de  Hulst  (p.  70)  for  details. 

In  this  case  the  cross -sections  n  duce  to  the  Rayleigh  approximation 


*1  1  XG 
2tr  3  X 


-  1 
+  2 


(9) 


25.  Kerr,  D.  E.  (Editor)  (1951)  Propagation  of  short  radio  waves.  Radiation 
Laboratory  Series  Vol.  13,  McGraw  Hill,  New  York. 
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or  Vail  Jo  Hulst  notation 


(!") 


(10) 


or  Van  do  Hulst  notation 

Q 


Q 


aH  irr' 


-  1 

4-  =  4x  Im  1  2 


in2  +  2 


(10') 


where  Im  denotes  "imaginary  part  of".  Defining  K  =  (m"  -  l)/(m“  -i  2),  Kqs.  (0)  and 
<  10)  become 


q  e  If  |K|2  - 

St)  J  v 


G 


(11) 


q  =  8rr2  Im  (-10 


(12) 


The  complex  index  of  refraction  m  is  defined 


m  =  n  -  ik 


<10 


where  n  is  the  index  of  refraction,  that  is,  the  ratio  of  the  speed  of  light  in  the 
particle  to  speed  of  light  in  air;  k  is  proportional  to  the  bulk-absorption  coefficient, 
Im  (-K)  in  terms  of  n  and  «  is 


Im  (-K)  =  - —2- -  (14) 

0  •}  2  2  2 
(n  +  k  )  +  4(n  -  k  )  +  4 

Thus  for  a  lossless  dielectric,  that  is,  one  in  which  k  -  0,  the  absorption  cross- 
section  is  zero.  In  this  case  the  extinction  (attenuation)  is  due  entirely  to  scatter¬ 
ing.  This  is  seen  by  Eqs.  (14),  (12)  and  (!')  where  Ini  (-K)  -  0.  For  a  lossy 


dulcet]  ic  a:;  wakT,  Hie  absorption  is  appre,  Table  .  In  f.n'1,  lor  x  l  tor  ivaliT, 
tin.-  scattering  d'lis:  ^eeliuii  is  negligible  in  comparison  wiili  tin-  aho  r-ption  i'mss- 

scct  ion  .uni  1.1  ,  a  Q  .  Dor  a  slightlv  lossjv  dielectric  as  ice,  Q  '  Q 

l'x,U  aj(  '  bi;  ‘‘i; 

(not  <  (inside  red  m  this  report). 

Tin-  die) ret l  it'  behavior  of  water  at  radio  frequencies  lias  been  icpo  le-d  m  the 
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literature  bj  many  authors:  Debye,"  Saxton,  “  Giant,  “  and  (.'ollie.  "  '1  he 

computation  prc>ce'dures  required  to  quantify  the  real  and  imaginary  terms  ol  the 
index  of  refraction  are  sueli  that  there  is  disagreement  among  the  values  calcu¬ 
lated.  In  order  to  avoid  these  computation  differences  (that  is,  round  off  errors), 
one  compulation  procedure  had  to  be  chosen  as  a  standard-  In  this  report  a  paper 
by  liay'^  has  been  chosen  as  tin  standard  because  Kay’s  research  covers  the  fre¬ 
quency  and  t  I'inpet  atui  e  range  of  interest  to  this  report,  and  Kay’s  mode!  ts  an 
empirical  model,  The  complex  refractive  indices  are  calculated  by  the  Subroutine 
Index  given  in  Appendix  11.  Typographical  errors  in  the  tables  uf  Kay’s  paper  nave 
been  corrected  in  this  subrout 'tie. 

Tin-  scattering  theory  discussed  above  is  scattering  from  a  single  particle. 

In  reality  fog,  clouds  jrj  tain  consist  of  not  one  size  particle,  but  many  sizes. 
Meteorologists  have  studied  these  dtopsi/e  distributions  us  depth  and  have  found 
them  to  vary  with  the  type  of  fog,  cloud,  ram,  geographical  location,  and  season 
of  the  year.  The  actual  dropsize  distributions  employed  its  this  report  are  given 
by  Silverman. 

Kaindcop  size  disti  ibutions  ate  the  most  complex  and  difficult  to  analyze. 

The  distributions  vary  temporally  and  specially  within  a  single  storm  due  to  the 
effects  of  wind  sorting  (updrafts),  evaporation,  accretion,  as  well  as  changes  in 
the  precipitation  process.  The  droplet  dial  ribut  ion  even  differs  in  different  parts 
of  a  single  rain  cloud,  both  vertically  and  horizontally. 

The  variation  of  dropsi/e  distributions  wil.lt  storm  type  was  shown  by 
I  'ui iu  arm  J  llis  results  showed  thunderstorms  had  a  much  wider  drop  size 
spectra  than  other  rain,  and  that  rain  showers  were  intermediate  between  con¬ 
tinuous  rain  and  thunderstorms.  These  results  are  borne  out  by  the  references 

.  „  .  32 

m  Carrier. 

As  stated  previously,  the  distinct  ion  between  clouds  and  rain  is  sometimes 
obscure,  especially  when  they  are  mixed.  Toy  and  clouds  are  assumed  to  be  com¬ 
posed  of  water  particles  which  are  less  than  iOO  pin  in  diameter.  For  particles 
this  Small,  air  currents  are  sufficient  to  overcome  gravity  and  the  particles  do  noi 
fall.  Clouds  are  generally  distinguished  from  fog  only  iu  their  distance  from  the 
earth.  The  fog  and  cloud  dropsize  distributions  are  given  in  many  original 

Due  to  the  number  of  references  to  be  included  as  footnotes  cm  this  page,  the 
reader  is  referred  to  the  list  of  references,  page  52. 
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research  papers.  Silverman  summarizes  both  fog  and  cloud  dropsize  distribu¬ 
tions;  his  table  is  reproduced  here  as  Table  2  with  an  additional  parameter,  liquid 
3 

water  content  ]\](g/m  ).  Note  that  the  mesoseale  parameter  (M)  liquid  water  con¬ 
tent  for  fog  and  clouds  are  the  same  order  of  magnitude,  whereas  the  microscale 
parameter  of  the  dropsize  distribution,  N  ,  <>  and  b  differ  by  orders  of  magnitude. 
On  the  microscale  level  fogs  and  clouds  are  distinguishable.  This  fact  is  impor- 
tanJ  for  infrared  calculations  of  transmittance  and  attenuation. 


Table  2.  Cloud  Models:  Drops i:.c  Distributions  and  Liquid  Water  Contents 


n(r) 

=  a 

r0,  exp[  -br 

Cloud  Type 

Model 

Number 

O’ 

b 

N  (cm 
o 

a 

M(g/'m3) 

Heavy  Fog  1 

1 

3 

0.3 

20 

0.  027 

0.  37 

Heavy  Fog  2 

2 

3 

0.375 

20 

0.  06592 

0.  19 

Moderate  Fog  1 

3 

G 

1.  5 

100 

2.37300 

0.  06 

Moderate  Fog  2 

4 

6 

3.  0 

200 

607.  5 

0,  02 

Cum  ulus 

5 

3 

0.  5 

250 

2.  604 

1. 00 

Altostratus 

6 

5 

1.11 

400 

6.  268 

0.41 

Stratoeumulus 

7 

5 

0.  8 

200 

0.  4rfGU 

0.  50 

Nimboscratus 

8 

1 

0.  333 

100 

11.  089 

0.  61 

Stratus 

9 

3 

0.  667 

2  50 

8.  247 

0.  42 

Stratus 

10 

2 

0.  6 

250 

27.  00 

0.  29 

Stratus- 

Stratocumulus 

13 

2 

0.  75 

to 

o 

52.  734 

0.  15 

Stratoeumulus 

12 

2 

0.  5 

150 

9.  375 

0.  30 

Nimbostratus 

13 

2 

0.425 

200 

7.  676 

0.  65 

Cumulus  - 

Cumulus  Congestus 

14 

2 

0.  328 

80 

1.4115 

0.  b? 
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■Silverman's  advcctive  and  radiative  fog  models  (termed  heavy  fog  and  mod¬ 
erate  fog  respectively  in  this  report)  are  in  the  opinion  of  the  authors,  the  most 
33 

representati ve  models,  Advective  fog  is  produced  by  the  transport  of  moist  air 
over  a  colder  surface  resulting  in  the  cooling  of  the  surface  layers  below  iiieir  dew 
points,  with  condensation  taking  place  in  the  form  of  fog.  Most  sea  fogs  are  ad- 
vective.  Both  the  size  range  of  particles  and  the  liquid  water  contents  at  e  large 
(see  Models  1  and  2,  Table  2).  Radiative  fog  is  produced  when  stagnant  moist  air 
is  in  contact  with  ground  that  has  become  progressively  cooler  during  the  night 
due  to  radiative  cooling.  The  cooling  from  below  produces  a  temperature  inversion 
in  the  layers  next  to  the  ground.  This  type  of  fog  has  both  a  small  size  range  of 
particles  and  a  small  liquid  water  content  (see  Models  3  and  4,  Table  2). 

5 

Cloud  droplet  size  distributions  are  more  difficult  to  model.  Mason, 

34  6 

Fletcher,  and  Borovikov  are  the  texts  consulted  by  the  authors.  Research 
4  32 

papers  by  Silverman,  Carrier  et  al,  and  the  references  therein  were  consulted. 
The  paper  by  Carrier  et  al  was  not  chosen  because  the  dropsize  distributions 
quoted  were  found  to  be  in  disagreement  with  original  references.  Table  3  is  our 
fitting  of  the  gamma-function  distribution  to  Carrier's  Table  1. 


Table  3 

Dropsize 

Parameters  for  Carrier's 

Data 

No 

Or 

A 

Ax  -  p 

cmode 

b 

a 

M  g/m3 

Stratus  I 

464 

8 

0.  85714 

2.  286 

10.  01 

0.  16 

Aitostratus 

490 

8 

1.0 

1.33 

4.  6 

0.  40 

Stratocumulus 

3  50 

3.  5 

1. 257 

1.0 

30.  09 

0.  24 

Nimbostratus 

330 

1.0 

2.  7 

0.  28 

26.  7 

1.4 

Fairweather  Cu 

300 

8 

0.  85714 

2.28G 

12,  68 

0.  1 

Stratus  11 

2  60 

3.  5 

1.  2667 

0.777 

7.  182 

0.  37 

Cumulus  Congestus 

207 

1.  5 

1.  914 

0.4285 

18.  71 

0.  43 

C  umulon  imbus 

72 

3 

1,4 

0.  6 

1.  555 

0.  16 

33.  Stewart,  D,A.  (1977)  Infrared  and  Submillimeter  Extinction  by  Fog.  TR-77-9 

U.S.  Army  Missile  Research  &  Development,  Redstone  Arsenal,  Alabama. 

34.  Fletcher,  N.H.  (1962)  The  Physics  of  Rain  Clouds,  Cambridge  University 

Press,  London. 
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Special  notice  should  be  given  to  Carrier's  liquid  water  contents.  In  a  private 
communication  with  Silverman,  the  large  discrepancies  between  the  liquid  water 
contents  calculated  using  a  gamma  function  and  those  quoted  m  the  literature  were 
discussed.  The  values  quoted  by  Silverman  are  typical  clouds,  based  on  data  from 
Borovikov  et  al^  and  other  sources.  The  variability  of  liquid  watei  "ontont  may  be 
ns  largo  as  a  factor  of  five  from  those  quoted  in  Table  2.  Thus,  it  may  be  seen 

how  difficult  it  is  to  model  clouds  with  representative  drops ize  spectra  and  liquid 

35 

water  content.  Research  has  shown  '  that  the  liquid  water  content  of  non-preeipi- 

3  3 

tating  clouds  have  values  from  0.  1  g/m'  -0.  5  g/m  whereas  precipitating  clouds 

3 

often  have  LWC  greater  than  1  g/m  .  Figures  9,  10  arid  11  are  dropsizc  distribu¬ 
tions  for  fog  and  clouds  (given  by  Table  2).  To  represent  the  dropsize  distribution 
in  Figures  9  to  11,  Diermendj  ian's  family  of  distribution  functions,  which  he  called 
the  modified  gamma  function,  are  used 

n(r)  =  a  ra  exp  (-br"1^  (pm  1  m  **)  (15) 

where  r  is  the  droplet  radius  in  units  of  pm  and  the  values  of  a,  y  are  related  to 
a  and  b 

b  =  —  -  (16V 


N-yb^/'r) 

r(ii  +  1  /y) 


(17) 


where 

rc  is  the  model  radius,  that  is,  radius  of  maximum  concentration 
r  is  the  gamma  function 

N  is  the  total  number  of  particles  per  unit  volume 


r„ 

c 

N  =  J  a  r<r  exp  (-br^)  dr 


(18) 


when  y  =  1,  =  0,  r g  =  »;  N  =  T(r)  thus  the  name  modified  r  distribution. 

85.  Blau,  H.H.,  Fowler,  M.G.,  Chang,  D.  T,  and  Ryan,  R.T.  ( 1972)  Cloud  micro¬ 
structure  studies,  ERT  P-375  Final  Report  NASA  Contract  NAS  5-21896. 
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RADIUS  (MICROMETERS) 


-3  -1 

Figure  9.  The  Dropsize  Concentration  (cm  pm  )  to 
tide  Drop  Radius  (pm)  for  Fog  Models  1-4  (Table  2) 


NUMBER  DENSITY  DN/DR 


flFGL.  CLOUD  MOOELS 


RADIUS  (MICROMETERS) 

Figure  10.  The  Drops  ize  Concentration  (cm  ^  pm  B  to  Par 
tide  Drop  Radius  (pm)  for  Cloud  Models  5-9  (Table  2) 


By  adjusting  the  four  constants  a,  b,  a  and  y  the  modified  gamma  function  can 
be  adjusted  to  fit  the  various  fog  and  cloud  dropsizc  distributions.  Tables  2  and  3 
have  7  set  equal  to  1  for  computational  simplicity.  This  follows  the  work  of 
Kingian  and  Mazin,  3<i  It  should  be  emphasized  again  that  the  fog  and  clouds  con¬ 
sidered  are  composed  of  water,  and  water  is  verj'  lossy,  that  is,  lias  a  relatively 
large  imaginary  component  of  index  of  refraction,  so  that  the  extinction  is  due  to 
absorption  only 


Q  »  Q  =>  Q 


EXT 


R 


O 

R 


The  meteorologists  usually  quantify  clouds  by  the  mass  of  condensed  water  per 
unit  oF  air 


00 

M  =  1  irp  J  n(r)  r3  dr  (g/m3)  (19) 

o 

where  p  is  the  density  of  water  in  units  of  g/'m3. 

Defining  the  volume  extinction  cross-section  for  all  the  particles  as 


'EXT 


00 

=  /  QEXTn(r)dr 
o 


(20) 


that  is,  the  integral  of  the  extinction  cross-section  over  the  distribution  of  drops, 

3 

Using  Eq.(19)and  Eq.  (12)  withp  set  equal  to  1  (g/m  ),  Eq.  (20)  becomes 


QV 

EXT 


Gn 

\ 


M  Im(-K) 


(21) 


The  result  is  that  in  the  Rayleigh  region  the  cloud  attenuation  is  due  to  the  total 
liquid  content  of  the  clouds.  It  should  be  emphasized  again  that  this  is  true  only 
because  of  the  frequencies  considered  in  this  report,  that  is,  the  size  parameter 
x  «  1.  The  size  distribution  (Eq.  (lb))  is  used  to  determine  the  liquid  water  con¬ 
tent  of  clouds  and  fog  via  Eq.  (19).  As  stated  previously,  the  meteorologists  quan¬ 
tify  clouds  by  liquid  water  content  (mesoscale)  not  dropsize  distribution 

36.  Khrgian,  A.K.  and  Mazin,  I.  P,  (1956)  Analysis  of  methods  of  characterizing 
raindrop  distribution  spectra,  Tr.  Tsents.  Aerolog.  Observ.  Moscow 
(translation)  17:3  6-46.  * 
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(microscale).  Thus  the  modelling  efforts  for  clouds  and  fog  is  much  simplified 
and  the  extinction  calculations  are  just  multiplications  (Eq.  (21)).  Function 
GAM  FOG  in  Appendix  D  accomplishes  those  calculations. 

The  calculation  of  attenuation  (extinction)  for  rain  requires  both  the  lull  Mie 
series  of  Eq.  (2)  and  the  dropsize  distribution  of  Eq.  (15).  The  adaptation  of 
Eqs.  (3)  through  (8)  to  a  computer  is  given  in  Deirmendjian  (p.  14).  This  pro¬ 
cedure  has  been  followed.  Deirmendjian  also  enumerates  the  useful  proportion  of 
the  modified  gamma  function,  specifically  Eqs.  (16)  to  (18)  above. 

If  the  Deirmendjian  distribution,  Eq,  (15),  were  to  be  used  for  rain  calcula¬ 
tions,  the  distribution  function  would  have  to  be  made  an  explicit  function  of  rain 
rate,  because  the  rain  rate  is  the  meteorologically  measured  paranwcer.  This 

procedure  of  relating  the  modified  gamma  to  rain  rate  (R)  is  shown  by  Dyer  in 
3  7 

Cole  et  al.  Another  procedure  would  be  to  use  Eq.  (25)  to  relate  R  and  the 
Deirmendjian  distribution.  In  this  report,  rather  than  basing  the  calculations  on 
purely  mathematical  (theoretical)  procedures,  an  empirical  dropsize  distribution 

function  will  be  used  which  relates  n(r)  to  rain  rate  (R)  explicitly. 

3  8 

Laws  and  Parsons'  (L-P)  provided  the  first  extensive  empirical  data  of  rain 

9 

rate  vs  dropsize  distribution.  Later,  Marshall-Palmer  (M-P)  proposed  a  drop- 
size  distribution  based  on  their  own  research  and  that  of  Laws  and  Parsons  which 
has  the  form 

n(D)  =  nQ  exp  [-bD]  (22) 

where 

n  =  8  X  lof  cm  1  m~3 
b°  =41R-°'21 
b  is  in  units  of  cm  1 
R  is  rain  rate  (mm  hr  *) 

D  is  diop  diameter 

Notice  the  original  definition  used  drop  diameter,  not  radius  and  n(D)  in  units  of 
-1  -3 

cm  m  .  The  units  of  dropsize  distribution  are  not  standardized  (see  Table  4). 
In  addition  weather  radar  specialists  often  use  N(r)  instead  of  n(r),  however  for 
consistency  we  will  use  n(r). 


37,  Cole,  A.E.,  Donaldson,  R.J.,  Dyer,  R,(  Kantor,  A„J.,  and  Skrivanek,  R.A, 

(1969)  Precipitation  and  Clouds,  AFCRL-69-0487, 

38.  Laws,  J.  O.  and  Parsons,  D.A.  (1943)  The  relation  of  raindrop  size  to 

intensity,  Trans.  Amer.  Geophysical  Union  24:452-460. 
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Tab] e  4,  Units  Utilized  foe  Scattering  by  Distributions 


Radar  Specialist 


n(r) 


or 


„  1 

rlD)  in  in 


-1 

cm 


Weather  Radar  Specialists 


n(r)  or  n(D)  in  m 


-3 


-1 

mm 


Optical  Specialists 


n(r)  or  n(D)  in  in  pm 


-■i 

n(r)  or  n(d)  in  cm  pm 


Meteorologists  and  weather  radar  specialists  often  use  dropsize  distribution 

-1-3  -1-3 

as  a  function  of  diameter,  and  n(D)  is  either  in  units  of  cm  ni  or  mm  m 

- 1  -3 

Optical  and  infrared  specialists  use  radius,  and  n(r)  is  in  units  of  pm  cm 
In  this  report  n(r)  will  be  used  witn  units  of  mm  1  m  3;  thus  Eq.  (22)  becomes 

n(r)  =  16000  exp  [-82R~0-21r]  mm'1  m'3  (23) 

The  changes  in  n  and  b  from  Eq.  (22)  to  (23)  are  the  result  of  changing  from  drop 
°  -1-1 

diameter  to  radius  and  units  of  cm  to  mm  .  Numerous  other  values  of  n  and 

o 

b  have  been  published  which  allow  the  M-P  distribution  to  fit  measured  data  more 

39  40 

accurately  in  varying  rain  (synoptic)  situations  (sec  Plank,  Joss  and  Table  5 
below).  Other  raindrop  size  distributions  have  been  employed  by  researchers 


n(r)  =  arc~4  exp  (-(r/b )A) 

-d  *41 

Best 

(24a) 

3 

n(r)  -  ar  exp  (-br) 

Khrgian  &  Mazin3(> 

(24b) 

n(r)  =  a  r  exp  (-br1/2) 

Deirmend’ian1  ^ 

(24c) 

39.  Plank.  V,  (1974)  Hydrometeor  Parameters  Determined  from  the  Radar  Data 

of  the  SAMS  Ram  Erosion  Program,  AFCRL-TR-74-0249TEnvironmcntal 
Research  Paper  No.  47?),  AD 

40.  Joss,  J.  and  Waldwogel,  A.  <19 50 >  Raindrop  size  distribution  and  sampling 

size  errors,  J,  Atmos.  Sci,  26:566-509. 

41.  Best,  A.C.  (1950)  The  size  distribution  of  raindrops.  Quart.  J.  Roy. 

Meteorol.  Soc.  76:i6-3G.  "  ’  ' 
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Thu  distributions  in  Eq.  (24)  contain  several  parameters  which  must  be  adjusted 
to  fit  empirical  data.  Unless  data  are  available  from  other  sources,  only  one 
parameter  can  be  fixed.  The  Marshall-Palmer  distribution  with  n^  fixed,  has 
only  one  parameter  b.  This  is  another  advantage  of  the  M-P  distribution  over 
other  distributions.  Even  though  many  exceptions  do  exist,  the  M-P  distribution 
appears  to  be  adequate  (see  Mason's  Appendix  B).  U  is  the  most  representative 
distribution  for  general  application,  and  like  all  the  distributions  given  above  it 
possesses  an  exponential  decrease  wit),  increasing  drop  size.  It  should  be  men¬ 
tioned  that  the  above  distributions  are  unimodal  even  though  there  is  experimental 
evidence  that  bimodal  distributions  exist  at  high  rain  rates.  A  point  specifically 
mentioned  by  many  authors  is  that  the  M-P  distribution  overweights  the  drops  of 
very  small  sires.  This  need  not  be  of  concern  because  the  range  of  integration 
employed  in  the  calculations,  starts  at  0.01a  nun  (see  Table  S,  section  b). 

Drops izr  distributions  are  difficult  to  measure.  What  is  measured  is  the 
rainfall  rate  (precipitation  rate)  H,  in  units  of  mm  hr  \  at  ground  level  with 
rain  gauges  or  at  elevated  levels  by  means  of  radar.  The  rain  rate  R  is  related 
to  the  dropsizc  distribution  through  the  integral  (see  Battan  p.  71) 


2 

R  =  ~  n  p  /  n(r) 
rl 


lv(r)  -  Vj  dr 


(25) 


where 

R  is  rain  rate  (mir,  hr  *) 

n(r)  dropsize  distribution  (mm  1  cm  “ ) 

r  dr  up  radius  mm 

v(r)  terminal  full  velocity  of  drops  in  (cm  sec  1) 

Vu  updraft  velocity  (cm  sec  *) 
p  is  density  of  particles  (liquid  water  density) 

Physically  rain  rate  is  a  flux;  its  units  are  g  cm  ~  sec  1  which  meteorologists 

convert  to  mm  hr-1.  In  calculations  it  is  usually  assumed  that  v  =  0  and  v(r)  is 

u 

given  empirically  by  Gunn  and  Kinzer,  Eq,  (25)  allows  an  internal  computer 
check  on  rainfall  rate  used  in  the  M-P  distribution  above.  Employing  Eq.  (23)  in 
Eq.  (2  5)  the  value  of  n  may  be  optimized  for  rain  rates. 

42.  Battan,  L.J.  (1973)  Radar  Observations  of  the  Atmosphere.  Chicago 

University  Press. 

43.  Gunn,  R.  and  Kinzer,  D.  D.  (1949)  The  terminal  velocity  of  fall  for  water 

droplets  in  stagnant  air,  J.  Meteorol.  6:243-248. 
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Table  5.  Characteristics  of  Precipitation  Derived  from 
Raindrop -Size  Distributions 


Source 

Location 

Type  of 
Precipitation 

No  (mm  hr  *) 

b  (mm  *) 

Marshall- 

Palmer 

Ottawa 

continuous 

8,  000 

4.1R-0’21 

Joss 

Switzerland 

drizzle 

30,  000 

5.7R"0,21 

cont  inuous 

7,  000 

4.  lR'0,  21 

thunderstorm 

1,  400 

3.  OR'0-  21 

Joss's  value  of  =  7000  gives  good  results  for  rain  rates  less  than  50  mm  hr 

M-P  continuous  and  Joss  drizzle  are  used  in  Figure  12  to  illustrate  the  LWC  of 

precipitation,  The  range  of  values  of  r  and  r„  for  various  rain  rales  as  given 
37  l  z 

by  Dyer,  in  Cole  et  al  are  indicated  in  Table  G,  section  a.  Section  b  values 
are  the  values  used  in  this  report. 


Table  6.  Range  of  Integration 


a)  Dyer 


Rain  Rate 

(mm  hr'1) 

(umi) 

i-2  (mm) 

R  <  1 

0.33 

1.  13 

1  <  R  £  5 

0.  03 

1.5 

R  >  25 

0.  75 

2.  3 

This  report 

R  (mrn  hr'1) 

Tj  (mm) 

r~2  (mm) 

.  25  £  R  £  150 

0.  015 

3.3 

As  seen  in  Table  6  for  the  rainfall  rates  considered  in  this  report,  the  range 
of  integration  over  radius  is  from  .015  mm  to  a  maximum  of  3.3  mm.  The  range 
of  rainfall,  rates  From  0.25  to  150  mm  hr  *  is  chosen  because  it  is  representative 


«r — - 
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of  rain  rail's  encountered  world  wide.  Ham  rates  as  high  as  300-400  mm  hr 
may  be  encountered  in  tropical  rains,  but  have  a  time  duration  of  only  about 
1  minute.  Haiti  rates  of  this  order  are  loo  specialized  to  be  considered  in  this 
report,  where  the  main  emphasis  is  on  world-wide  atmospheric  models. 


RAIN  RATE  (mm/ii) 


Figure  12.  Liquid  Water  Content  (g/rv>‘  )  vs  a  Function  of  Rain  Rate  (mm /hr) 
for  Drizzle  Employing  the  Joss  Model  for  Dropsize  Distribution  and  for  Rain 
Employing  the  M-P  Model.  The  relations  are  Joss  M  =  .  HR-  88;  M-P 
M  =  . 072R- 88 


5.  MODELS 


This  section  will  review  typical  mid-latitude  fog,  cloud  and  rain  models.  The 
general  concepts  for  such  models  are  given.  It  should  be  re-emphasized  that  the 
theory  of  scattering  by  hydrometeors  is  well  known,  and  that  the  meteorological 
statistics  of  these  scattering  particles  are  not  so  well  known.  This  fact  impacts 
the  systems  designer  of  millimeter  and  submillimeter  systems.  The  completely 
deterministic  propagation  system  modelling  of  the  microwave  region  has  to  be 
foresaken  for  the  statistical  propagation  models  of  the  millimeter  and  submilli¬ 
meter  regions  forced  on  the  system  by  the  physical  properties  of  the  atmosphere. 
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In  somi'  systems  the  luxury  of  "w  .si  ease"  .analysis  for  propagation  is  not  pc  iri- 
hie  because  lor  the  worst  ease  the  atmosphere  may  be  opaque  for  a.l  practical 
purposes. 

It  is  a  simple  procedure  to  identify  the  physical  and  meteorological  pt  rameters 
required  for  calculation  of  attenuation  and  transmission.  They  are 

n(r)  -  dropsizc  distribution  of  hydrometeors  which  is  related  to  liquid 
water  content  M  and  tlie  l  ar.ifall  rate  It 

X  -  wavelength  of  incident  electromagnetic  wave 
m  -  index  of  refraction  of  the  drops 
T  -  thermoinetrie  temperature  of  the  drop's 
h  -  thickness  of  the  fog,  cloud,  rain  system 

In  addition  it  is  necessary  to  know  the  frequency  of  occurrence  of  clouds,  and  when 
modelling  a  cloud/rain  situation,  it  is  necessary  to  define  a  consistent  model  of 
cloud  liquid  water  content  and  rain. 

The  parameters  n(r),  X,  m,  and  T  have  been  reviewed  previously.  The  latter 
three  parameters  of  the  above  list,  that  is,  thickness  h,  frequency  of  occurrence 
of  clouds,  and  a  consistent  water  content  model,  will  be  expanded  on  in  this  section, 

5.1  Cloud  Models 

Clouds  are  classified  as  low,  middle  or  high  according  to  the  following  rough 
44 

estimates  by  Berry  et  at, 

low  clouds  0-2000  meters 

middle  clouds  2000-0500  meters 

high  clouds  >  6500  meters 

In  addition  to  these  classifications  vertical  development  of  clouds  should  be 
included.  Tabic  7  relates  the  general  classification  of  clouds  to  specific  cloud  type 
nomenclature. 

In  this  report,  only  water  clouds  are  considered.  These  clouds  belong  to  the 
families  of  low  and  middle  cloud  classification.  High  clouds  (cirrus)  are  usually 

r.  r 

composed  of  ice  crystals.  References  for  cloud  types  are  Mason,  "  Borovikov, 

44  32  21  45  46 

Berry  et  al.  Carrier,  I.uke,  Diem,  Weickman  and  .Aufra  Kampe, 

47  48  49 

Durbin,  Gates  and  Shaw,  and  Squires  and  Twomcy. 

Cloud  water  content  is  related  to  cloud  droplet  spectra  by  Eq.  (19).  Table  2 

as  originally  proposed  by  Silverman  is  for  "typical"  clouds.  These  "typical1 

clouds  are  not  average  values.  For  example.  Cloud  Model  5,  cumulus,  has  a 

Due  to  the  number  of  references  to  be  included  as  footnotes  on  tills  page,  the 
reader  is  referred  to  the  list  of  references,  page  53. 
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liquid  water  content  of  1  g.m  (  fable  2),  This  type  of  cloud  may  have  typical 
Values  ol  liquid  water  from  0.  a  to  1  g  in’  with  values  as  high  as  ■}  g  in'  depend¬ 
ing  on  geographical  lo>  alien  (lor  example.  New  ling  land  vs.  Florida),  Othvi  types 
ol  cumulus  clouds  may  have  liquid  water  eoiilenis  as  high  as  8-10  g  in'*.  Liquid 
water  eeiileiit  of  clouds  is  very  important  this  study  Localise  the  cloud  droplets 
are  Un\leigh  seatlerers  (see  Figures  G  to  tl,  Mie  vs  Kaylcigh,  expanded  scale,  and 
temperature  dependenee,  respectively  and  Figure  Id,  which  in  addition  to  showing 
the  temperature  dependence  of  the  r  Tractive  index  explain  the  temperature  n'oss- 
ove.'  of  Figure  8). 

The  frequency  of  occurrence  of  clouds  and  cloud  types  is  beyond  the  scope  of 
the  present  paper.  Systems  designers  and  reseats  hers  should  consult  the  follow - 

fit)  f)l  r  ■)  I) 

ing  relerenees:  Solomon,  Cult  man,  Quaylc  el  ai,  Visher,  and  Jlorovikov. 


Table  7.  Clouds  Classified  According  to  Heights 


Cloud  Heights 

Low  Clouds 

Stratus 

St 

0-2000  meters 

St  ratocunmluK 

St' 

Nimbostratuy 

Ns 

Middle  Clouds 

(C 

1 

o 

o 

A  Host  rat  us 

As 

1 

Alto-umulus 

Ac 

(not  considered  in  Table  2) 

Vertical  Development  Clouds 

Cumulus 

Cu 

700-8000  meters 

Cumulonimbus 

Cb 

700-20,  000  meters 

(not  considered  in  Table  2) 

50.  Solomon,  I.  (19(13)  Estimated  Frequencies  of  Specified  Cloud  Amounts  within 

Specified  Ranges  of  Altitudes,  AWS  TH~No.  1(177’  "  ~ 

51.  Gutlman,  N.J3.  (1971)  Study  of  Worldwide  Occurrence  of  Kog.  Thundershower! 

Supercooled  Low  Clouds  and  Ft  equeney  Temperatures,  N avail'  50-1C-G0, 
distributed  by’NV/SED,  Asheville,  n7~G1 

52.  Quayle,  11.  C..,  Meserve,  J.R1.  and  Crutcher,  71,  L.  ( i 9GG)  Climatological 

Mean  Probability  of  Penetrable  optical  Path,  distributed  by  NOA^b" 
Asheville,  N.C. 

53.  Visher,  S.  S.  (10 (50)  Climatic  Atlas  of  the  United  States,  Harvard  University 

Press,  Cambridge,  Mass. 
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Figure  13.  The  Ileal  and  Imaginary  Components  of  the  Complex  Dielectric  Con¬ 
stant  (Calculated  from  Ray)  for  0°C,  2Q'-’C  and  40°C 


The  two  cloud  types  altocumulus  and  cumulonimbus  of  Table  7  are  not  con¬ 
sidered  in  Table  2  because  too  few  drop  size  measurements  have  been  made  in 
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these  clouds.  Gates  and  Shaw  consider  nimbostratus  and  fairweaLher  cumulus 
to  bracket  altocumulus  clouds. 

In  the  millimeter  and  submillimeter  frequency  range,  since  clouds  arc 
Rayleigli  seatterers,  altocumulus  can  be  modelled  by  assuming  a  liquid  water 
content  between  that  of  the  nimbostratus  -  and  fairweather  cumulus  (that  is, 

-^•0.5  gm  . 

Table  H  relates  the  cloud  models  of  Table  3  to  cloud  thickness.  The  range  of 
values  of  thickness  may  be  large  for  these  models.  For  example,  advectivc  fog 
ran  be  from  30  meters  thick  to  300  meters  thick  depending  on  meteorological  con¬ 
ditions  and  geographical  location. 
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Table  8.  Typical  Cloud  Model  Thickness 


Silverman's  Model 
Number 

Height  (m) 

Bottom  Top 

1 

0 

1  50 

2 

0 

150 

3 

0 

75 

4 

0 

75 

5 

660 

2700 

(i 

2400 

2000 

7 

6G0 

1320 

a 

160 

1000 

9 

160 

660 

10 

330 

1000 

1 1 

660 

2000 

12 

660 

2000 

13 

160 

660 

14 

660 

2700 

660 

3400 

5.2  Rain  Models 

Precipitation  as  it  affects  propagation  requires  a  knowledge  of  the  variation  of 
rainfall  rate  and  dropsize  distribution  with  altitude.  Such  data  have  to  be  related 
to  the  surface  measurements  of  rainfall  rate.  The  occurrence  of  various  intensi¬ 
ties  of  rainfall  cannot  be  obtained  directly  from  climatological  records.  Usually 
the  climatological  data  are  available  for  average  monthly,  seasonal  and/or  annual 
totals  of  precipitation.  The  quantification  of  rain  rate  (per  hour)  is  subject  to 
much  uncertainty  because  rain  gauge  data  are  recorded  every  6  hours.  Through¬ 
out  Europe  and  the  United  States  clock-hourly  precipitation  data  are  available  at 
selected  weather  stations.  Clock-hourly  rate  is  the  total  precipitation  on  the  hour 
every  hour.  The  instantaneous  rates  of  precipitation  may  vary  considerably  within 
the  clock-hourly  rate.  For  example,  a  0.  2  mm/hr  precipitation  reported  during 
a  clock-hour  might  have  accumulated  in  a  30  minute  period  at  a  rate  of  0,4  mm/hr 
or  in  a  10  minute  period  at  a  rate  of  1.2  mm/hr.  Tipping  bucket  rain  gauges  may 
provide  instantaneous  rain  rate;  these  data  are  available  at  specific  weather  sta¬ 
tions.  In  addition  to  climotological  data,  weather  radar  data  are  available  to 
describe  the  temporal  and  spatial  range  of  rain  rates  at  areas  in  Europe  and  the 
United  States. 
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Two  standard  precipitation  ii. tensity  classifications  arc 


1.  Standard  Weather  Chart  Code 


Code 


Hate  (in/hr) 


Ha  in 


•  • 
.% 


(intermittent  light  rain) 
0-01  in/hr 
0.  n-0.  3  in/hr 
>0.3  in /hr 


continuous  light  rain 
continuous  moderate  rain 
continuous  heavy  rain 
intermittent  moderate  rain 


2. 


Intermittent  rain  is  not  quantified  in  inches  per  hour. 
The  National  Weather  Service  Weather  Radar  Code 


Code 

0 

1 

2 

3 

4 

5 
<; 


Rate  (mm/'hr) 
0 

0-3.  8 
3. 8-20 
20-40 
40-70 
70-152 
>  152 


Rain 

none 

light 

moderate 
strong 
very  strung 
intense 
extreme 


The  modelling  of  a  physically  realizable  cloud/rain  system  with  altitude  is 
extremely  difficult.  The  model  has  to  be  consistent  with  respect  to  liquid  water 
content  (that  is,  the  cloud  associated  with  a  specific  rainfall  rate  must  have  a 
l'quid  water  content  the  minimum  value  of  which  is  of  the  same  magnitude  as  the 
liquid  water  content  of  the  precipitation).  For  example,  the  cloud  associated  with 
a  5  mm/hr  surface  rainlall  rate  must  have  0.  25  g/m  cf  liquid  water.  Water 
storage  in  clouds  is  another  factor  which  has  to  be  considered  in  modelling.  There 
are  three  stages  to  a  rain/cloud  system,  the  buildup  stage,  the  fully  developed 
stage  and  the  decaying  stage.  These  stages  differ  in  total  liquid  water  content. 

We  have  attempted  to  model  the  fully  developed  stage  only,  that  is,  the  models 
have  maximum  liquid  water  contents.  The  models  assume  maximum  updrafts  for 
maximum  water  storage  above  the  cloud  base.  In  other  words  these  are  "worst 
case"  conditions.  For  applications  to  "typical"  situations,  one  can  safely  assume 
the  total  water  content  of  the  rain/cloud  system  is  one-half  that  given  m  the  models 
(see  Table  10).  The  difficulty  in  modelling  rain/cloud  systems  is  that  the  researcher 
does  not  know  which  stage  of  development  the  system  is  in  unless  measurements 
are  made  in  real  time.  Also,  the  liquid  water  content  can  vary  from  0.  5  to  5 
g/m'  depending  on  the  stage  of  development,  season,  geographical  location,  etc. 
Thus,  it  is  obvious  that  only  tile  fully  developed  rain/cloud  stage  can  be  successfully 
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modelled.  When  applying  the  models,  one  has  to  keep  in  mind  that  they  are 
"worst  ease"  models.  "Typical"  models  require  one -half  the  values  quoted  in 
Table  10.  In  addition  to  the  cloud  liquid  water  content,  the  precipitation  liquid 
water  content  beneath  the  eloud  base  has  to  be  modelled.  Table  0  relates  liquid 
water  content  M  (g/m3)  as  a  function  of  rain  rate  II  (mm/hr) 


Table  'J.  Relationship  Between 
Liquid  Water  Content  and  Rain 
Rate 


M  -  0.  072R '  88 
M  =  0.  071V  83 
M  -  0. 052R- 97 


b4 

Table  9  is  Atlas's  data  taken  from  Sissenwine  and  it  applies  to  both  cloud  and 
precipitation.  Other  relationships  have  been  used  in  which  the  exponent  of  R  is 
set  equal  to  1. ,  lor  example,  M  =  0.05R. 

In  addition  to  the  liquid  water  of  clouds  and  rain,  the  water  vapor  in  clouds 
should  be  assumed  to  be  at  saturation  (100  percent  humidity). 

Drizzle  and  orographic  rain  cannot  be  modelled  as  described  above.  The 
modelling  of  a  drizz.le  is  questionable  because  a  drizzle  is  comprised  of  a  large 
amount  of  small  droplets  the  diameters  of  which  are  of  the  order  of  100  pm.  In 
this  size  range  of  particles  ;t  is  difficult  to  separate  cloud  from  precipitation. 

The  drizzle  for  a  stratus  type  of  cloud  system  has  a  rain  rate  <1.  25  mm/hr  or  a 
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liquid  water  content  of  -0.  10  g/m  and  the  cloud  associated  with  the  drizzle  may 
have  IAVC  as  high  a  1  g/m3.  This  associated  stratus  cloud  may  be  in  contact  with 
the  ground  or  as  high  as  500  m  above  the  ground  depending  on  the  meteorological 
conditions.  These  considerations  have  been  taken  into  account  in  the  drizzle  model. 
Drizzle  Model  1  is  a  typical  drizzle  (see  Figure  14),  Tt  is  possible  to  model  drizzle 
as  a  precipitation  rate  R  up  to  the  cloud  base,  and  then  a  LWC  in  the  cloud  which 

Q  O 

is  the  sum  of  two  equations  M  -  0.05K  (cloud)  and  M  =  0.0721V  °  (precipitation). 
These  same  equations  have  been  used  in  Models  II,  III  and  IV  (see  Figures  IS 
through  17).  Table  10  relates  model  altitude  to  liquid  water  content. 


54.  Sissenwine,  N.  (1972)  Extremes  of  Hydrometeors  at  Altitude  for  MIL-STD- 
210B,  AFCRL-72-0369,  AD 
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Figure  17.  Model  IV:  Summer  Cumulus,  15  mm /hr 


Table"  10.  Rain  Models 


Unit 

Distance 

liquid 

Water  Content 

(km) 

Model  Name 

Base 

Top 

Fully  Devol 

jpccJ  "Typical" 

I. 

Drizzle 

0.  0 

0.  5 

0.  55 

0.  27 

<1 .  25  mm/hr 

0.  5 

1.0 

1.  10 

0.  5f] 

1.0 

1.  5 

0.  25 

0,  12 

II. 

Steady  Rain 

5  mni/hr 

0.  o 

0.  5 

0.  55 

0.  57 

0.  5 

1.0 

1.  65 

0.  83 

3<R<8  mm/hr 

1.  0 

2.  0 

2.  20 

1.  10 

2.  0 

3.  0 

1.10 

0.  55 

3. 0 

4.  0 

0.  30 

0.  15 

nr. 

Steady  Rain 

0.  0 

1.0 

0.  66 

0.  66 

12.5  mm/hr 

1. 0 

2.  0 

3. 87 

1.  93 

7<R<20 

2.  0 

3.0 

5.  16 

2.  58 

3.  0 

4.  0 

2.  58 

1.  29 

4.0 

5.  0 

0.  77 

0.  38 

IV. 

Sumtn  or 

0.  0 

2.  0 

0.  80 

0.  80 

Cumulus 

2.0 

3.  0 

4.  65 

2.32 

15  mm /hr 

3.  0 

4.  0 

6.  20 

3.  10 

4.  0 

5.  0 

3.  10 

1.  55 

5.  0 

6.0 

0.  93 

0.  46 

The  Models  II,  III,  and  IV  as  given  in  Table  10  have  been  generated  from  data 
'H  as  hfi  ^7 

of  Valley  and  the  research  of  Wexler  and  Atlas,  East,  Kessler,  Dangille 

and  Gunn,  ,lfi  and  Hamilton.  J-'>  fn  order  to  derive  general  concepts  (rules -of-thumb) 

to  model  rain,  the  research  of  the  above  authors  has  been  compared  with  the  M-T 

formula  and  t lie  cloud  liquid  water  content  formulae.  The  factors  required  to 

multiply  the  M-P  formula  and  the  cloud  LWC  formula  are  given  in  the  following 

steps: 


55.  Wexler,  K.  and  Atlas,  D.  (1955)  Moisture  supply  and  growth  of  stratiform 

precipitation,  J.  Meteor.  15:531-539. 

56.  East,  T.W.R.  (1957)  An  inherent  precipitation  mechanism  in  cumulus  clouds. 

Quart.  J.  R.  Met.  Soc.  S3jGl-69. 

57.  Kesslet,  R.  ( 1 9 G 1 }  Kinimatical  relations  between  wind  and  precipitation 

distributions  II,  J,  Meteor,  13:510-. 

58.  Langiile,  R.C.  and  Gunn  K.  L  S.  (1948)  Quantitative  analysis  of  vertical 

Structure  in  precipitation,  J.  Meteor.  5:310-. 

59.  Hamilton,  P.  M.  (I960)  Vertical  profiles  of  total  precipitations  in  shower 

situations.  Quart.  J.  Roy.  Meteor.  Sc-c.  92:346-3  62. 
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1.  The  liquid  water  content  of  the  ulou<l  is  expressed  by  the  formula 


M  (g/m3)  0.  OSH 

N. B.  for  I!  <  8  mm/hr  the  cloud  (d  <  100  pm)  may  touch  the  earth's  surface,  as 
in  Model  II.  For  higher  rain  rates  (H)  this  is  not  true  (see  Model  111);  at  these 
higher  rain  rates  the  larger  drops  scavenge  the  smaller  droplets. 

2.  At  the  cloud  base  the  sum  of  the  LV/C  of  the  rain  and  the  cloud  is  multi¬ 
plied  by  n  factor  of  3  for  a  unit  altitude.  For  Model  I  the  unit  altitudes  are  all 

O. 5  km  For  Model  II  the  first  unit  altitude  is  0.5  km  and  subsequent  unit  altitudes 
are  1  km.  For  Model  III  the  unit  altitude  is  1.0  km.  For  Model  IV,  a  summer 
cumulus  rain,  the  cloud  base  is  raised  to  2  km  and  the  cloud  top  to  6  km. 

3,  At  an  altitude  of  one  unit  distance  above  the  cloud  base  the  LWC  of  the 
rain  plus  cloud  is  multiplied  by  4,  the  next  unit  distance  the  LWC  is  multiplied  by 
2  and  the  last  unit  distance  the  LWC  is  multiplied  by  0,  (>. 

Figure  18  is  an  example  of  tha  attenuation  over  a  1  kilometer  horizontal  path 
through  clear,  cloudy  and  rainy  atmosphere. 

It  should  be  emphasized  that  the  above  models  are  to  be  used  only  as  examples 
of  typical  rain/cloud  models.  The  computer  program  given  in  Appendix  D  can  be 
used  to  model  any  real  meteorological  system.  It  is  the  responsibility  of  the 
researcher  to  determine  what  is  a  realistic  model  for  the  geographical  location 
being  considered. 
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NEPERS/KM 


FREQUENCY  (GHZ) 


Figure  18.  Attenuation  vs  Frequency  for  a  Clear  Standard  Atmosphere,  a  Cloud 
of  0.  18  g/m3  LWC,  and  a  2.  5  mm/hr  Rain  Rate  for  a  1  Kilometer  Horizontal 
Path 
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Appendix  A 


Definitions 


a)  Transmittance  and  attenuation 

As  stated  in  the  text,  the  AFGL  line-by-line  absorption  calculation  code 
HIT RAN,  was  generated  by  infrared  specialists.  These  specialists  use  trans¬ 
mittance  to  characterize  atmospheric  propagation,  whereas,  microwave  special¬ 
ists  use  attenuation  to  characterize  the  same  atmospheric  propagation.  The 
attenuation  (a)  is  related  to  transmittance  (r)  by  the  formula 

a  =  1  -  r 

The  transmittance  at  a  single  frequency  for  a  path  of  length  x  is 

_  emergent  power  _ 
incident  power  '  o 

For  optical  on  IK  propagation,  the  rate  of  attenuation  is  generally  given  as  an 
attenuation  coefficient,  (j,  which  is  expressed  as  the  fractional  reduction  in  the 
transmitted  energy  per  unit  length  along  the  direction  of  propagation 


so  the  units  arc  reciprocal  length  there  km'1) 
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b)  Decibels  and  Nepers 

I'm1  microwave  and  radio  propagation  the  attenuation  is  generally  given  in 
terms  of  a  decibel  (dB)  or  neper  (less  commonly)  loss.  Decibels  are  defined  by: 


number  of  decibels  -  10  log 


10  l  I 


=  10  log10  T 


A  1  neper  attenuation  means  a  reduction  to  l/e  of  the  original  value  or  number  of 
nepers  =  In  (I/I  )  -  (number  of  decibel$)/4.  343 


/3(km  *)  -  |3(neper/km)  -  |3(db/km)/4. 343 


Optical  depths  are  dimensionless.  It  should  be  kept  in  mind  that  both  nepers 
(logarithms  to  base  e)  or  decibels  (logarithms  to  base  10)  are  d imensionless  units 
and  the  two  are  related  by. 


INp  =  10  l.og10  e  dB  -  4.343  dB 


It  should  be  noted  that  neper  loss  is  frequently  applied  to  the  voltage  or  amplitude 
of  the  electromagnetic  signal  on  v  ave  and  not  to  the  transmitted  energy  on  power, 
as  we  use  above.  This  leads  to  an  extra  factor  of  2  relating  nepers  to  decibels. 

1  neper  reduction  in  amplitude  =  8.  686  decibel  reduction  in  intensity 
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FUNCTION  GAMFOG 


Appendix  B 

Subroutine  AERSOL 


calculates  attenuation  due  to  fog/cloud  liquid  water 
content  (g/m^)  by  the  equation 

7f  =  8.  1888  l  M  Im  [- 

X  \  m2  +  2 
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FUNCTION  GMRAIN 


where  A  =  wavelength  in  cm 

3 

M  =  liquid  water  content  (g/m  ) 

m  -  complex  index  of  refraction  calculated  in 
subroutine  INDEX 

-  calculates  extinction  due  to  rain  by  the  full  Mie  theory, 
as  a  function  of  wavelength  temperature,  and  drop-size. 

-  drop-size  distribution  is  Marshall-Palmer 


n(r)  =  N  EXP(-br) 
o 


-  DATA  WVLTAB  is  wavelength  table  used  in  ATTAB 

-  DATA  ATTAB  is  attenuation  data  table  in  dB/kni 

for  temperature  20°C  Table  B1 
(abridged  version) 

-  DATA  TMPTAB  is  interpolation  data  table  for  tem¬ 

peratures  in  °K  Table  B2 


-  DATA  TLMDA  is  interpolation  data  table  for  wave¬ 
length  in  cm 


-  DATA  RATTAB  is  rain  rate  table  in  mm/hr 


Subroutine  INDEX 

Applied  Optics  vol  11  No.  8  Aug.  1972  pp  1838-1844 
Typographic  errors  in  original  paper  have  been  corrected  as  follows: 
pg  1837  A  2:  8000  p  ,  =  2(t  I  273)/5  +  72 

should  be 

A  &  300  cm  ,  c  =  2T/5  -i  88 

pg  1840  -  Table  11  A-  column  sixth  from  top  0.08  should  be  0,009 

1  -4  -S 

-  equation  12  d.  10  should  be  10 

Table  111  -w  column  fifth  from  top  688.24  should  be  588.24 
°i 

-  Equation  nvr)  =  n^  (7,0  -  A)  +  n^  (A  -  6A) 

should  be  n(r)  -  nH  (7.  0  -  A)  +  n1-1  (A  -  6) 
r  r 
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Appendix  C 

Computer  Input  Instructions 


Programs  required:  FASCOD1,  HIRACl,  LBLF4,  SCANFN1,  FCODMRG  and  the 
subroutine  package  AERSOL. 

To  utilize  the  Aerosol  Subroutine  in  calculation  of  Attenuation  in  the  Microwave 
region  (1-lO'JO  GHz)  the  input  control  cards  must  be  as  follows: 

*5.  1  CARD  1:  XIDU),  1=1,  7) 

FORMAT  (10A10) 

AN  IDENTIFICATION  CHOSEN  BY  THE  USER  TO  BE  PRINTED  AT  THE 
START  OF  THE  OUTPUT  AND  AS  A  HEADER  FOR  ANY  PLOTTED  FILE. 

5.2  CARD  2:  IHIRAC,  JLBLF4,  ICNTNM,  IAERSL,  IEMIT,  I3CAN,  li-LOT, 

MPTS,  NPTS 
FORMAT  (715,  35X,  315) 

IHIRAC  =  1  CALL  HIRACl 

VAN  VLECK-W EIS5KOPF  APPROXIMATION  TO  LINE  SHAPE 
USED  IN  MICROWAVE  REGION 
IHIRAC  =  2  CALL  HIRACL 

LORENTZ  APPROXIMATION  TO  LINE  SHAPE 
IRIRAC  =  3  CALL  HIRACD 

DOPPLER  APPROXIMATION  TO  LINE  SHAPE 


The  description  of  input  control  cards  1-7  is  brief  and  accurate  at  the  present 
time.  For  a  detailed  and  up-to-date  description  see  FASCODE1  Report  by 
Clough  et  al.  8 
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ILBLF4(0,  1)  Flag  for  I.Bl.Fi  (LINE  BY  UNIS  FUNCTION  -1) 

1  -  Yes  (extends  calculation  to  258  o’  ) 

ICNTNM  =  0,  1  Flag  Tor  continuum  (CONTNM) 

1  =  Yes 

IAEKSL  -  0  No  AHRSOL  Calculations 

*  1  AERSOL  Calculations  desired 

(MUST  FOLLOW  INPUT  OF  AERSOL  DA  S A  AS  EXPLAINED 
IN  SECTION  5.  9) 

1EMIT  =  0  No  Emission  -  Only  Attenuation  (Absorption  Coefficient) 

=  1  Emission  and  Transmission  calculations  performed 

ISCAN  (0,  1)  Flag  for  SCANFN  1  =  Yes 
IPLOT  (0,  1)  Flag  for  PLOT  1  =  Yes 

MPTvJ  -  0  Number  of  points  to  be  printed  from  panel 

NPT£  =  1-20 

Controls  number  of  output  data  points  to  be  printed  from  merge 
routines. 

5.3  CARD  3:  VI,  V2  (cni'M 

FORMAT  (2E10.3) 

Beginning  and  ending  wave  number  values  for  tile  calculations, 

5.4  CARD  4:  TBOUND  (°K) 

FORMAT  (E10.3) 

Boundary  temperature  (for  ICNTR1-1  CASE  only) 

5.5  CARDS:  ANGLE  (degrees) 

FORMAT  (E10.3) 

Secant  angle  as  measured  from  local  zenith.  Calculation  performed  for 
Angle  of  0°-00°  and  ANGLE  of  120°-180° 

0°-60°  sets  control  (ICNTRL-3)  for  looking  up  case 
120°-180°  sets  control  (ICNTRL=1)  for  looking  down  case 
otherwise  secant  =  0,  and  ICNTRL=0. 

IF  ANGLE  NOT  IN  THESE  RANGES  -  SECANT  AND  iCNTRL  WILL  BE 
READ  FROM  LAYER  CARD  (see  card  7) 

5.  G  CARD  8:  NLAYRS 

FORMAT  (15) 

NUMBER  OF  LAYERS  INPUT 

5.7  CARD  7:  PAVE,  TAVE,  SECNTK,  1CNTRK,  (FILHD(I),  1=  1,  4) 

FORMAT  (3F10.  4,  15,  4A10) 

PAVE:  Average  pressure  (MB)  for  the  layer 

TAVE):  Average  temperature  (°K)  for  the  layer 

SECNTK:  A  non  zero  value  for  SECNTK  will  override  the  value  calculated 

from  Angle. 


j 
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ICNTRK  -  II-  *■  0  It  replaces  value  d  ICNT11L 

=  I,  Looking  from  HI  toward  112  for  111  >  112 

-  2,  Looking  from  III  toward  112  for  III  --  112 

=  3,  Looking  from  Ill  toward  112  for  Ill  <  112 

(PILIID(I),  I  I,  4):  Descriptive  text  defining  beginning  and  ending  of  layer 
in  (KM)  (Col.  30-75) 

5.  8  CARD  8:  (1VK(M),  M - 1,  7) 

FORMAT  (7E10.  3) 

WK:  Array  containing  the  molecular  densities  for  the  layer  in 

(Molecules /cm^) 

for:  1I20,  C02,  Of ,  N.,0,  CO,  CI14,  C2 

5.9  After  each  input  layer  of  CARD  7  and  CARD  8  and  cte.  The  following 
procedure  must  be  followed  in  order  to  utilize  the  Aersol  calculations. 

5.9.1  CARD  9:  NLYRS 

FORMAT  (15,  21-TO.  3) 

NLYRS:  Musi  be  1  less  than  number  of  eld,  fog,  and  or  rain 
cards  read  in. 

NLYRS  1  Blank  or  0:  signifies  no  cloud,  fog,  and  or  rain 
between  that  particular  layer. 

5.9.2  CARDS  10-NN:  Where  NN  is  11-29 

The  maximum  number  of  input  CLOUD,  1-00,  and  or  RAIN  data  eards 
between  each  main  program  layer  is  20.  If  more  a"e  needed  than  change 
the  dimensions  of  RRA,  llDIl',  T2  and  C2  in  subroutine  AERSOL  as  desired 
11,  P,  T,  HI1,  CLOUD,  RAIN 
FORMAT  (F5.  1,  5F10.3) 

H:  Height  of  Layer  in  KM 

P:  Pressure  for  that  Height  m  MB 

T:  Temperature  in  (°K) 

RH:  Relative  Humidity  values  of  0  —  1.  0 

CLOUD:  Liquid  water  content  of  cloud  or  fog  in  (g/m^) 

RAIN;  Rainfall  rate  in  (MM/HR) 

5.  10  Example  of  card  squencing  for  a  specific  case. 

Problem:  Calculate  Absorption  Coefficient  for  three  2  KM  layers  with 
clouds  and  rain  in  the  first  4  KM  for  0-40  (cm  wavenumbers. 
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(*  REPRESENT'!  A  A  PACK  ON  TI1E  CARD) 


CARD  1: 
CARD  2: 
CARD  3: 
CARD  4: 
CARD  5; 
CARD  6; 

CARD  7: 
CARD  8: 

CARD  ft: 
CARD  10: 

CARD  11: 
1st 

Layer  CARD  12: 

CARD  13: 
CARD  14: 


EASCODE  MICRO  WAVE  REGION  TEST 
**  !>?1  1 

***0.  0***  ***-**40.  0 
BLANK 
*****()_  () 

****3 

**903.  25**281.  7624*******************0.  0 KM  to  2.  OKM 
2.  fl25E4  22*l.  517134  21*1.  3 55E4  12*1.  287E+18 
*3.  448E+17*7,  357+18*9.  G32E+23*3.  590E+24 
***4 

**0. 0****1013. 3 *****291. 2*******1. 0*****0, 000 
********2.  00 

**0.  5*****954.  6*****288.  9*******1,  0*****C.  015 

********!_  91 

**J. .  0*****898.  8*****28G.  7*******1. 0*****Q,  04  5 

******** 1^  79 

**1 . 5*** **345. G  *****284, 3  *******1 . 0*****C. 078 
******** 1 ^ 

**2.  0*****795.  0*****282.  2***-****l.  0*****0,  106 
********1# 43 


2nd 

Layer 


CARD  15:  **705.  0873**208.  3 198****<t»************2,  0KM  to  4.  0KM 

CARD  1G:  *1.  24.'E+22*i.  245E+2t*l.  252E+17H.  056E+18 

*2.  828E+17*6.034E18*7.  901E+23*2.  945E+24 

CARD  17:  «***4 

CARD  18:  **2.  0*****  795.  0*****282.  2*******1.  0*****0.  106 

y*******! , 43 

CARD  19:  **2.  5*****746.  {>*****279.  8*******1.  OWW^O.  134 

**#**#**!_ 25 

CARD  20;  **3. 0*****701.  2*****277.  7**#****i,  0 ^>*^**0.  158 

********1. 06 

CARD  21:  **3.  5*****G57.  8*****275.  3*******!.  O*****0.  176 

********0, 87 

CARD  22:  **4.  0*****61 6.  G*****273.  2****#**1 . o*****0.  160 

********0, 70 


I  CARD  23: 

3rd  CARD  24: 

Layer 

I  CARD  25: 


**543.  7893**255.  8787*******************4.  0KM  to  60KM 
*4,  531E+2DM.  010E+21*1.  142E+X7+8.  5S8E+17 
*2.  295E+1 7*4.  89GE+18.  G.  41  lE-t  23*2,  390E+24 
ulatik  (no  rain  and  or  cloud  this  layer) 


LAST  CARD:  #.  0001*****’!I**0.  o******40.  0 


Appendix  D 


Listing  of  Program 


SUBROUTINE  A^RSOL 


rit/TU 


0PT»1 


FTN  <u  7+M71 


i 

5 


Iff 


15 


20 


25 


30 


30 


*.0 


45 


50 


55 


C  * 
C  * 

c  * 
c  • 
c  • 

c  * 
c  * 

G  * 
C  * 
C  * 
C  * 
C  • 
c  * 
c  * 
c  * 
c  * 
c  * 


ccc 

ccc 

ccr 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 

ccc 


SUBKOUn^C  *FPSCLURAD> 


PROGRAM  rtEPSOL 

CALCU.fltrc  ATTENUATION  IN  THE  HICROVAVE, MILLING  TER 
ANT  Sn-HIU-lMETER  RFGtONS 

1  ►'AR  1*379 


*  a/thqfs 

* 

*l  L  EON®  9  f  W.  A9REU  A • F .G. L .  617-061-4775 

* 

*  VUCENT  j.  FALL  ONE  A  •  F  ,  G  *  L  *  aUTOWON  475 -4775 

9 

*  #♦****•*******•♦*•'****•*********♦****#***•*»¥♦#.♦»**♦*#¥  *  *  «*«#*•»' 

CCHMON  /a  95rRP/ViO, V2D,0VC,WPT0,0 (205 ) 

COMMON  /FI  l  HOF/  AID  (ID  )  ,  SECANT  ,  PA  V£  »  T  AVE  ,  HMCLIO(  2  0  >  ,WK  <20)  ,MN2  ,0V 
IVlf  V?,TnO  UFO,  f  SCO  ID  (17)  ,NHOL  ,NLAYCR 
CIHENSION  TPF (7051  ,RRA<20) ,HOIF  t20> 

DIMENSION  T?(?0>  rC2(?0> 

NLYRS  4  J  S 1  PE  1  LESS  THAN  NUMBER  OF  CLO  AND  OR  RAIN  CARCS  INPUT 
pftheev  each  l a  TER  pead  in  from  MAIN  PROGRAM. 

blank  r«R.n  RF«p  IN  -  SIGNIFIES  NO  RAIN  AND  or  CLO  DATA  acThEFN 
THAT  PMTirULAR  LAYER. 

FROM  H-  SJrOF'O  LAYER  ON  IT  IS  NFCESSARY  to  repeat  tfe  previous 
TCP  LF/“ L  INPUT. 

EG.  I(PUT  AFTEP  main  program  later  OF  0-2KN  is 
7  =NLYRS 

0.0  KM  ETC.  ♦  . 

1.0  KM  ETC-  .  . 

2.0  KM  ETC.  .  . 

THEN  I*  «>UT  A  FT  £  R  MA  IN  PROGRAM  LATER  OF  2-4  KH  IS 
2  sNLYRS 

2.0  KM  r  TC.  ,  . 

3.0  KM  eTC.  ,  . 

4.0  KM  ETC.  .  . 

MAXIMUM  NUMBER  OF  INPUT  CLO  ANO  OR  RAIN  LAYERS  BETWEEN  E ACM 
MAIN  PROGRAM  LAYER  IS  70. 

REA  0  93  5  *  N  L  YRV 
IFCNLY7S.  NF.O)  GO  TO  9fl 

return 
90  CONTINUE 
TP*  C  .  0 
CP*  C. 0 

READ  yiO^HjPtTjRMpCLOUDtRAlN 
PRINT  ?19»HfP,T,RHpCLOUO,RAIN 
OLDT*T 

olon*h 
FIRSTHM 
OLOC*ClOUO 
OLOR*R*  IN 
PRINT  318 
OO  110  L*1 *NLYR$ 

REAP  910»MfPfT iRHfCLOUO t RAIN 


68 


SUBROUTINE  *ERSPL 


7i4/7i  0PT*1 


FTN  4I.6HU0 


0 


PRINT  919»H,°|T,RH,  CLOUD  ,  RAIN 
^•'i  r  (L)  sM-rtop 

NO  R*A  <L)  =  <0.  TP+RAIN)  /  2. 

T2<1 )  =  (OLf>T*T)/?, 
r?(LI-tPL)r+CI  OUD ) / ?• 
cs ( l  >  =n?(. )  *Hnr  tu 

TP  =  1 F  *T2 IL  )  *  C  2 ( L ) 

55  cP=rp*r?<i> 

PRIM  9  24,  (’L0H,HrTL  (L)  ,  C  ?  (L  )  ,  5  RA  (  L> 

CLP  1  - T 
CL  PP=H 
OLOC=ClCUt 

7  D  0LPP=WTN 

110  FONT  1  ivlii  r 
F  INfiL  H=M 

TRIM  91*>  FTRSTH,  FINAL  H 
f  IPCTH-  j 

75  Tftv=TfVril 

PET  A«l.  r  .6951  *1.1  v) 

CALL  SFCONTL^O) 

CO  ?l»l  Js’jNOTO 
VJ=  vio+rt"*  M  ( J-l)  *OVP 
90  C J- VJ*?9.q 7925 

rasr LOr  0.0 
DO  12 Z  L=1 , NLYRS 

fAMLLO=r,A-f(LO+GANFOG<GJ,T?<L>  ,  C?(L>  > 

95  G9RANTrr.1*nUfH0lF(U*GNR«lN<GjfT2<U  ,  R  R  A  (  L  )  ) 

iro  rcNiiNiif 

CJ-GMPA  NT/  til  *+& ANCLO  /4 ,  343 

c°f  <j)-rj 

tXPV.I-f  YP(  -VJ*PETS1 

90  pAPF  0C=  VJ*  <l.-EXPVJ>/<  l.*EXPVJ) 

IMJ.FQ.7I  0  ’=PJ/RADFAC 
IF  (  j.t.0,  *)  DT^DJ/RAPFAC 
IKJ  )  =0<J>*  rj/RAOFAC 
201  CONTINUE 

95  CALL  SECONrtxn 

TIKE  =X1-X0 

CM  )  =C<  1J  ♦r?-(03-D?> 

PRINT  9i%  Tint 
PRINT  9i7,Vl0,V?0 

100  PRINT  920,  <0<J> ,J*1 ,NP TOf 

301  CONTINJE 

905  FORKATI 15, ?F1 r ,3) 

910  FORHATl  F5.  1  f  5  F  1 0  •  31 

915  FORHATLIOX,*  TINE**  ,F1  0.  5/) 

105  91G  F0RKATC1M1  ,10X,*HGT=*,  f  5*1#*  TO  %F5.1//> 

91  7  FOPKATC  15X  ,*V10=*,F  10.  3,  *V20=*,F  10,  V/5X|*0  (J)  ,J=l,KPTO  RAIN  ♦  CL 

inur*//> 

918  FORKATdMl  ,5X,*  INPUT  DATA  OVG  TEMP  AND  AVG  CLOUD  L YR5*/ / ) 

919  FORMAT  C//1TX,F5* l  ?5  FlO • 37/) 

110  920  FORMATE  5X,  10C1  2.5) 

924  FOR  HAT  M  0  <  *  *  OL0H=¥  ,F5  .  1  f*  M**,F5.1»*  T2**,F10.3i*  C2=*,F1Q.3, 

1*  RkA=*,Fl  0.  V) 

RETURN 

ENO 


1 


THIS  PAGE  IS  HEST  HUiUU-J-i 
KHIOM  GUP*  Fl&SttSKMJ)  10  HD.C 


TVI3  A. i 'Cl  QA  KL*R 

1  AVK\A*^v-r-:  “ 


G9 


i 


v 

* 


T 


■~S'0?K  COTTnE  t'NfiPV 


0PT-=1 


TTN 


SUnprdT^e  lNlUX<MVl  rfTC,  KEY  ,  RE  AL  ,A  I  m  r: ) 


*  *  *  >  *  *  *  * 


~Z~* 

c  *  * 

“c  '•  *  wAvrur^CTH  is  i.j  rt>T imp Teas.  ^f*pf:9atu*£  is  in  deg.  c«  . . “*  *' 

C  «  * 

T“  *  •  kfy  ;s  "5'fTTo  t  in  SUfrouTinc  T^mfo  g  ~ 

r  •  * 

‘  r.  ♦  ■»  RFfli  i 


1h?  RFAl  PART  OF  THl  ftEFkACTm  TmCFa. 


L  *  *  M  M  A  C  TS  r  ME  IMiCINA.I  PARt  Cf'TML  f?£  F  fiftC  TI  \r  l  TNCEX  IT  I'l 


"*-$7T;rfi  ms  n  r-sru  “  'T*‘ma/r 


c  *  * 

C  *  •  A  $  L  R  If  S  of  CHttKS  Are  VSflE  '4'k-0  wARrJIMJS  I-TVi?1!, 

c  *  * 

C  *  *  r<AV  APr\.TfO  0°  T  ICS  VOL  ll,N0.8,Alj£>  ?*,  FT*  ld't-lP<»<i 

C  *  * 

-  r  '  •—  rvj 6 £*£ ' fCTfS"  u.v/£  °3 E7T'h  PDF  T”c  <Tm  Y  “ S*  f  0 TS ~\~£Tpl  ^ 

C  #  M 
C  *  * 

«*****>•  *»»¥•.*•*»*»  »«».».  *  + 

c  i  -  l  • «: 


V  * 


75 


Tr  <  wa  v.  TTTTTTcoi)  wp|TE1o*1i  ’  ~ 

imr.i  r,-?M  wpitf(6,?> 

1  fo?pat  r///f  33x,  *  aTTemfT  i  nG  Tc  €vAuja?E  FCR  A  kA  VElEnGTH  lFSS  ThAM 
j-tjf  Mirpr.N** //) 

FCrrATT ///^AY  ;>AYTr»iFTlNG  T  C  tVAUJATF  F0C  *  Tr*PF&ATURt  LFSS  TF  Ah 
7  -’f  .  rEGRTES  CLNTIGRAr -:*,//) 

77, l  L  Dc  F yp  (WAVL,TF,KF  y  ^YAlT*  1  *7^' 

»■  *  T  A  nL  F  3  WAT6D  PG»  1?*«C 

<.  lr<wAv*  .FT,,D>i>'  r,7  . 

F  i  F  <  VlA  VI  i  f'T^i  i)  1  1  ,  F 

l  ??  "=fi:r  (wAvC‘,V.li3ACo  ,1  .4  ,  lo3«.^.2  .'2^  »?45005, , 

_ ?.  ,161  *29  ,4371°, 7,  27Ff5t  7>  _ 

kC  A  I*  ft:  AL*  (WL^L-.r3<*>  /  .  QG6  *P  l  *  <  .  1  -W  A  V  U  )  /  .  06  F 

- GO^rrl  1 - 

f  IF  (HA  VI  .CT..OOOF)  e,lC 

■  r«Tnn  to  vrrmyiwri-6Tr.T^  2  rorrca  - 

4’5<i?l  3,  , 161 . 29  ,4371^. 7, 27G61  ,?) 

- mi=Rr  AL*ftc«r*crr-?  -r.  >  * .  h 

7F  (  HA  VI  «  GT  •  .  GC  Q  7 )  tl,  9 


rcwsrnOTo  - — " 

^ ~irr=T"0P (TOVCVirJ'^TT  ,"3^2  .27  i?.19F3F^  B<»7T»-^ ^TS.^r 


r  c: 


n-»  r  tj  r 


,  24  f«7 


:  a; , 


,7C1 A  ?DSi 


R"l= PiVk  i*€Tft»'25.  \  *.'Ci&Ci*Sx'p<  t .  0iOO2v^  wAvLV”  •'2^> 
RFAL=R1  *  ( .  COOT-HAVU/  .  3  0  0 1  ♦«£  At  *  <  «A  /L-.0PC6)/.000l 
ro  TO  1  1 


50 


10  RrT  A  Lk  DCF  (HA  VL,  1. 709  07  ,3  752.27,99 .914  F  «0<tf  15 .196  3E+  0  4, 1639,  , 

- STTCtt F 3 79^*4 G. ?7,F6fi .TV^B' V 74 E TTE'f T *Tf  10  77B1 5rfD 5T 

RFAL=Rf  ALsRFAL*  <Tr-2t»*)»,0001»eXP((  •  OOCoaF’KAVL)  »y.2  5) 

*  -■»  TABUE  2  “TOTPTT^i  f^TTg  ” 

11  IF(K7’M  •  C  f  •  •  3 )  GC  TC  5/ 

- rrrovn  SFmrTnTrrc - 

1?  AIMAG=A  IflAG^  4B  CWAVL  25 , 300,  ,•  47f  3.  U  «B(MAVL,.39f  17.  f  .45  f  1.3>  ♦ 

- 7tot  w  vi— ^  r,  visyi  m - - 


GO  TO  5  7 


•nu 


s  ?k" 


;  «-> 


,v  -  v  v  1 

xu  - 


70 


■sirairDtrrrwTfroET 


— 74  770 - OFTVT) 


■o 


13  IF<H*VI  ,GE..  006?)  14,15 

- it;  "7siinrG^riT3)iT;*)iBTiiT»cTr<iTVii7vrr35iirn'>n^3Wvn-m  ,777774571  .37  *~ 

63  Be  fl<  Ha  Vl  ,  .  25,  300.  ,  .*1  ,2.) 

_ — — - - fro  TO  "67 — - - - - - 

15  IF<  NOVI  .GF..0017)  16,17 

r^TTHTo :-T T1- sru  am  wj  vl  ,.  .’9 ,  f7. 65, it 3  n  eBnrm: rnsi  j&r,  ,722  ,i.*m 

9 18 (  we VI  ,  .25,  30  0.  4,?.) 

T5 - tnrr<rT7  - - - 

17  IF<  RAVI  ,G£.. 1)0  061)  10,19 

- — ’ - n  mre s=? rwi G+mTfrmmT? , 6 ny. o 47,  to imra vl ,  .79 , 77. ,vi65,7.T7 > 

1  a  Fj  <  Wfl  VL  ,  .41,  62.,.  22,1.0) 

jTfl— jtj  t7  —  — 

70  19  IFOiavi  .GF..  000495)  20,  21 

~  73  JT8eG  =  SimG*7PTV*Viyr0i  ,4793, 73571.  J  FABIHAVL  ,717,571 ,  .03  9,  ?.l — 

GO  TO  CT 

- 21  TFT  TraLjr.T.r.TD  D  073  75  77777 — ~~ - - - - - - 

2?  eiHeGs,IBeGF68(HeVL,.27,2.9  7,.0  4,2.»  +  eB(HaVL,.lil,  4.9  5,.0  6,l.) 

T9 - - - - - 'nrryT - — — - - 

23  A  IMAGED  IHAG*AO(WAVL  ,.£7,2,97,.  025,2.)+«B«HaVL,.01,4.  95,.06,l.) 

- 57  roNTTWE - - - 

RETURN 

- ENP - - - - - 


Sl'ERTCTJTrNF  TJFBYF 


7xmr — opt  =y~ 


TDPuWFI 


1  SUBROUTINE  0EBTE(HAVL»TC,KEY,RE,8I) 

- T7TT737T. -  “ 

IF<  KEY. KE.0)  1  ,2 

- - - rTFT S=E.  771  ?7+.T)?IEI,7!,TTC-.7T)131 1  93*TC^TC - 

5  *LPHe=-16.B1297T*. 0609265 

t  a  tjftbw  TrmrT5*nrp  1 2  err.  9  *m  —  1  1 

SIG  =  12,  5664E»08 

- F5=73 .<■  4» ri  .-7074F79*TTC -7571T. 0 0371 1 9*TTt-2r.T  »»2-rO00 333 37 8 »  ~ — 

1  1TC-25.  )  **3) 

id - — — — r<rTcn»  - — - 

EFI N=3*  168 

573'?ior7uonz3'TC*TT>  •  IT35,7,T(' 77733 
SIG=1  .?6*EXP«'12500  ./(T*1.9869>) 

- TAIJ^3.TT93mF-0TTF)fPa777l)7TrT*1.3*63VT  ~  “ 

15  FS=3.1C8  «.15*lC“TCl2.  5»TO20C. 

3Ti7T't'iro>T4Vl  — —  - 

f  2=  1.  57  0  8*  ALPHA 

- 277  »T7»3TT»TiyTC3TrjT_‘3TNTC7Tyrrr*T7r*TTT;riTP'HaT'T - - - 

risFFIMt-  l£S-EFIN)»»l.  *<Cl»*Cl,-ai.PHA)  *S  I N  (C2  )  »  )  /  OEM 
79  TH  *ET. .  97ffl  0  .  R  e  VL7G  F  .T33. 0  Fi=37*Tr/37*«»;  —  - 

F2=  lES-EFIMl'Cl'1*  <1  .-ALPHA  >*COS  (C  2)  7DEM»SIG*hA  VL/18.  0496E*10 

- . - - T1FS5QRT7TF3  *  SORT  fET^n >E2  *F7TJ777>  ‘ 

Af  =  -E2/<2.‘'RE> 

- - - frriT^j - 

25  fnp 


71 


TOTRrrramjoT’ 


7  **774 - (TPl  =T“ 


TTN  G  ,fc  WTfll - 1 


1 


S 


FUNCTION  OOP(HAVL,A  ,CEN1  ,  P  ,  C  ,  CE  N2  ,B  ,  £  ,  CE  NI,  F  ,  Gt 

W,W  vt  '  ■  -  ■ 

«?=vv 

*  *2^V2 - - - 

W2*CEN?**2-V2 

"F3*TEN7**2-V5  ' 

COPsGOTl <A+B*Hl/(Hi *H1»C*V2> *0*  HE/<  H2*HJ*E*  V2  M-F  *N3/ (M3*  H3 +G*  V2  )  > 

Trrtretr  — 1  • 

fno 


~TVm - 5PTST 


~FTTr~!i.'^g6D  — — ~ — U<i 


1  '  F U Nf  T 10 W  AB(WAVL  ,A,CEN,B,C) 

- ipi-»«F  xri^nwsi  mot-i « mrnnr;  *  wivt7CEN77B>T»*CT 

RETURN 

pjjjj  — 


~r^7Th - DFT=r 


THIT.TOn  “  TJt 


i 

FUNCTION  G«NFOG(FPOG,T,RHO) 

. “C 

c 

COMPUTES  ATTENUATION  OF  CONDENSED  HATER  IN  CLOUCS  OR  FOG 

5 

X, 

c 

FRO-  s  FREQUENCY  <  GH  ?) 

c 

c 

T75E?  Hfi  VETIUTOEP  “TIHVERSE  TTP~J 

T  =  1 1  MPt  Hfi  J  UKF  (LitOHtfcS  KELVIN) 

---  c 

c 

“PTTO  =  CONDENSED  JTATFR  5TATTJ R  CONVENT  m/CUVTlrH HE  TERT 

CTNOEY  =  COMPLEX  DIELECTRIC  CONSTANT  fl  FPGH  INDEX 

in 

t 

COMPLEX  ClNaFX,riNOX 

GftMFOG-0, 

15 

<  tTTJRN 

2  f  REQ  =  FPQG/ 29.979 

V.AVL  =  2°  .979/FRQG 

JT3.2 

C*Ll  irOEMWAVl  , TC, *E Y, REAL,AIHAK) 

70 

C IN  Dt  X=c  rpiTX  (R t  A  L  *  ATM  DX) 

GAMFaG=e,18  6d*FREQ^fiH0»AIMAGt  -  <C INDEX** 2-1)  /  (C.TNOEX**  ?+2l) 

'HETiJR’K  " 
FNn 


7 


7  4/  74  PTToi 


F  !  M"ltTb  '<t  b1] 


"0 


FUNCTION  GMRAIN<FREO,T,PATE> 


cowutes  attenuation  of  ccnoenseo  hater  in  form  of  rain 


FREO  =  FREQUENCY  <GH7> 

=  TEMPERATURE  I  Dk’SRFES"TriVTH I ' 


T 

RAT! 


RREC  IPITA  TI  CN  RATF  (HH/lfO 


TIT 


C  UTILIZING  MARSMAU.-FAi.MER  SIZE  DISTRIBUTION  HITH  RAYS  INDEX 

~c - bt-refehtito — - 

c 


- PWtT»TT?“W»Fi;TtiCTh  TlTOl'E  rOTrWWtFUCTRS  TJStO  IN  T  ABL”E  ATT  HIT - 

- tnat  i  rrs'  ivrrciFTnnrT^TiTnnrnTnnrTFHPTst  TijRET'TirTinrrt'rvTTr 


!  TOSTEtrHGTH  “Td  CR  ' 


IEITTAB  IS  RITV1  RTiT  TABUmnWHlT 


TA8LC  MIAg  FOR  REPRESENTATIYE  RAINS 


25 


nr 


“TTTKCT  TMTRPOOTTSL  SCHERFTTRTTrrN  ~9Y~ 
F.T,  FLORANCF  O.N.R.  PSSAOENA  CA. 


35 


nr 


45 


SB 


J  TTAITDR  T57E.  5  J 

OIHENSTON  *(3>,Y<3>,ATTN<3> ,RATES«3I 

ETH  ER  ETTTRr  TRRTABTSTTTTHD  A  IBT.T  ACTTTO  .TTAETTST - 

CATAHVL  TAB/  .03,  .033,  .0375,  .(143,  .  05,  .  OE,  .  V  T  5  ,  .  1 ,  .  1  5  ,  .  Z,  ,  ZS  ,  .  3  .  .  5  , 

1  ■  b  ,  1  a  ,7  »»  3  ,TG  .,5,  *5, 5  ,T)  •  ,1,  •  5  *  /  ■  ,  0.  »  9  •  ,111  a 


,15a/ 

pATA  RATTAB  /.  Z5  ,  t .  Z5  ,2 . 5 , 5  . ,  1Z  .  5  ,  Z  i>  .  ,  50  .  ,  1 00  .  ,  150  .  / 

— tnrrm  no*/,  03  , .  17.t;t.Z5,;:.:ntt>.7 
CATA  -  T*P7A  °/27  3  .  ,  26  3,  ,Z9  3.  ,  303.  ,  3 13  .  / 

d ATATisrnrB it ,tt  ,  i=tvz  7i  /,  urn  *  »btt 7.1104 , . 39  j  ,  •  9 0 r,  79 t  2  , .ns 9 , 

1.874,.  7  7  3,.  656,.  5  3  9,.  434  ,.179,.  06  34  , 

-  2.0’Si..  E85r-7r77m^  ,  TTTTr^rTErzr-T  ,VE7*r-  3  774T4  E  -  37 - 

3.363  E-f  ,.309t-3,.2*Ze-3,.18iE-3,.  14  5E  -  3  ,  .  6  3 1 E- 4/ 

- rTTTTATTVHTT  ,?1  ,  T=  r»77T7T7?  T7f«  YJ7Z  Tf,  7TTB7ZTRT,T  ,RT77TTT,-  ^ 

1Z,51,Z.  A1,Z.ZZ  ,1.99,1  .T4,.B19,.3Z'*,  >  _ _ 

3air7E-Z,.i3ZE-Z,  «4b  Yt  “3 ,  ,  Y41E-3,.  58BE-7,,245E-3/ 

- OTf  A  I  AT  T  A B  i  r7TT7TST7?5T  ZTITT 73 .95737(70  ,S  .65  ,TTZT;T.  77,3.  84, - 

1  3.90,  .3.  0.3,3.  63,3.34 ,3  .0  1,  l.TZ,  .783, 

— rrv!T7a  jor,  .onr;  ,ti  r<r,  .not  -tttvzte-tytt  74 1 -27 - 

3  .30  4E-E  ,  .Z5ZE-Z,  .18  Zi-Z,  .130E-Z,  .10  9E-Z,  .45  4E-3/ 

- CAT  A  v  AT  1  HP  074  )T7*T7  ?  75  /  5  .  J57T3.  4175741  75.5  6  75. 65, 57  76,  5.U9  , - 

16.01,6.  00, 5,  53, 5.  49  ,5.06, 3.  .9, 1.60, 

— JTF7.T 35 771750,  .  C ?!C  ,  7C 13  J'.TTl IT  -Z ,~ TSfr^TJ - 

3. 60  EE -Z , .495E~Z,.35iC-Z,.262t-Z,.Z04E-2,.BZ9E-S/ 

- DAT  A! AT  I  ABU, El ,  1=1 , Z ZT/9 . 35_, 9i 44TVTT6 ,T. 70 i’BTnFjnrTTSiTOTJO’i'" 

110.  59, 10. 80,10.69, 10.  S3, 9.  81, 7. 13  ,3. 94, 


55 


tunction  gw  rain  "7  u^t*i  mnrrswfi 


2  2.7  0,.6  9*..245  ,.  0  9ZO,.0  399,.02B2  ,.0  20  9, 

- - - V.ll  fcTV.  6TTST'.'*  57€  -5T.  t  i  1  f-?r."i7JT-T,Tr J5TT/  ~~  ' 

fit  OAT  A  (ATTAB(I  ,61  ,1=1  ,2  7)  71*  .  27 , 14 .40 , 1  4  .5  0 ,14. 79, 15. 0  3, 1  5.  33,  15.71, 

i'167ITi)  (THT? ,  16  .70  ,  lil  3*,  lTTJ  1,1?.  36, 7.  Si  , 

25.38,1. 52, .591,. 236,. 100, .0686, .0488, 
S.0S62..C279,.0li0,.0127,.S47E-2,. 34*E-2/  ~ 

_ OAIMATTABIJ  ,7)  ,7  =  1  ,27)  721. 78, 21. 97, 22. 24, 22, 54, 2Z. 90, 23. 36, 23. 9o, 


■S'? - 

- 1747  6*, A  5761,25  7**7*577*  rf-T.14 , 2  0  .1T7T37*77 

210.  3  7,?.  23,  1.3*,.  60  2,  .2  65,.  179,  .124  , 

3  *  tfS  ^2  yTTTS^l ,  *TTT9  9  ,  •  0Zc7  ,  ■  01 9? ,  •  F61E-Z  7 

CAT  MAT  TAB  U,B  ),  1=1,27)/ 33. 22,  33. 49 , 53-0 6 , 34 » 32, 

34.  .5,35*59,36.40, 

70 

T37755 ,79. 18,59. *4,  397  96, 39. 50, 34. 54, 24. *TJ 

219.40 ,f .  26,3.0  9,1. 4®, .70  6, .4  66, .338, 

“3 .239,  .1  72, . (ronrTTTSTr, .  04lb,.0 17*7 

CATS (AT TAB  <1 ,9 J , 1*1 ,2  7)742.48,  4  7,82,43,  3  0,43.84, 

44.51,45.53,46,47, 

147.93yC0.16,5l«  ill  ,51*54,51  71?  2,  .94(34.46, 

2  27.  59,10.0  6,4.  86,2.  4  9,1 .24  ,.872  ,.613, 

TV 

3.435 ,  .TTT|  . irrrnffj,.  0677,  .0198/ 

CAT  A  FACTOR  / 1,0  ,1.0, 1.0, 1.0, 1.0, 0,9  9, 0.99, 1.0,1. 

01 ,1.02, 

ii.rs, 1.01,1.1, i7*.f. o  .*? .  i .  jz ,  i  ro  ,17  o ,  *.99 , 

21.5  5,  1.  25, 1.0, 0.81,0.  £5, 1.72,1.  29, 1.0, 0.79,0  . 64, 

60 

31*0,1  of  ,1*0, 1  *  0  y  1  *  0  yl  .0  yl.flyl.O  ,1.0,1  .Dly 

41. 01, 1.  01, 1. 0,0. 9  9,  C.  98 ,0.95, 0.96, 1,0, 1.05, 1,10, 

5172*717  14, 1.0,0. *6 ,11 . 72 ,  777  371. 70  ,1 . 070 .797  0.66, 

61*0)  l»fj  1 «  0  »  1 •  0  y 1 • 0  y  1  fO  ,1.0  |1 t  C  yl.0yl«01y 

7 i  .TJ ?  , 1 »  0  i  1 1  •  di  D«4?i  Co  97  t  0 .  ?fc  1  ff*^7  y  1  •  0  y  1»  0 4  ,1.07  y 

81.  04,  1.  03, 1.  Q,  0.  95,  C.  88 1.30,1*1,  0.79,  0.63, 

85 

91.  S,  1.*, T7H,T7B,T7T,i  .0  ,1 .0  ,T.  9  .0  ,1  .*1 , 

11.02,1.01,1.0,0.98,0.  97,0.99,0.99,1.0,1.02,1.04, 

?''.91,C.96,i.0yl*0i,J. 01*1. 75 ,1.31, 1*0,0. 76, 0.62, 

31.0,1.0, 1.0, 1.0. 1.0,1. 0,1. 0,1. 0,1. 0,1. 01, 

91 

41.03,  1. 01, 1*  0  y  0  •  S  8  y  0  •  57,1*01 yl  *  0*  1*0,  1*0, 1*01, 

5  0.88, 0.  95, 1.0,  1,04,1.  06, 1.72, 1.31, 1.0, 0.78,  0.62/ 

- TATA  RBTLIH  7.05/ 

C  GlVt  ZERO  ATYN  IF  RATE  FALLS  BELOW  LIMIT 

TP <  RATP •  GliffATl. IN J  CO  to  iz 

GHRAIN  -  0. 

-7* 

RFTUF  N 

12  WVLT:»  *  29. 979/f REQ 

IF<WVl.iH.GT.HVLTAB(l>>  GO  TO  14 

ISO 

- ICCW* 

JH*7=2 

TOT* 

CCC  THIS  00  LOOP  IS  2  LFSS  THAN  NO.  OF  WVLTAO  INPUT 

14  ra  is  T=*r?5 

tF<NVLTH.LT.<.S*(HVLTAB(I»®MVLTAB«M>)l>  GC  TO 

16 

nr? 

■  1 - ir^oNTirot 

CCC  SET  ILCH  EQUAL  TO  1  LFSS  THAN  pn  RAN 

GO  TO  1  8 

no 

IE  HlOW  *  i-z 

16  CONTINUE 

00  190  12777 

IF  (RATE,  LT,f,5*(RATTAB(I)+RATTA0CI+l))))GC  TO 

195 

- - -  niircoHTTNut 

KHX  N*6 

i 


75 


Tuvcriofnnnwiv 


7 1,  m  onv 


rmirB«fFir 


115 

CO  to  <  98 

i *n  iw-j-i  - - - — - 

198  CONTINUE 

Co  ? U  o  —  1 1 JM  A?  _  - - - 

tj  =  no«  *  j 

MJ)  =  -  A  L OG1 H  VL  TRRTTJ7  i  - - - - - 

20  CONTINUE 

?  =  -01  OG(HL  LTH) 

125 

oinr'TiT-  ■  ■  - - - - - 

*J=  KH1’'  +  K 

WA1  fcS  <*  »  sRATTAB  <KJI  - -  - — - 

CO  2<i  J  =  1 .  JM  AX 

Tj  -  II  OH  +  J  ~  - - ■ — 

71 J)  •-f  LCG(STTA0)I  J  ,K.!i  ) 

1  JU 

24  CONT  I  i.t’E  - - - - - 

0TTN(K»=;Eyp<-0XTK(X,y,7tJP0X»  ) 

_ C_ 

2b  CONTI Wl't  "  '  __  -  -  - * - 

APPl  Y  TFHPERATURE  CORRECTION 

135 

IFFT.LT .THPTAB (I) )  CO  TO  37 

Jl  CON  TINt'E  '  ~  - - - - - - - 

TLOW  =  4 

k  id  ;i  ....  ■■■  - - _ — _ 

33  ILOW  =.  1-1 

00  PI  J  =  2.6 

IrlHwL'r.LT*  TL  NOfl  ( J)  J  CC  10  43  " - ~ - — - 

41  CONTINUE 

145 

ccc  JLOW  IV  Z  LESS  THAM  Uu  Hi  1  - - * - 

JL0W=4 

CO  TO  «*  5 - - — — - — 

43  JL  0  W  -  J  -£ 

4^  CONTJNI'f  - - - - - - - 

PO  50  7=1,2 

c 

uo  49  J=  1 ,  ?  - — - - — - 

INTFHFOLATE  IN  TEMPERATURE 

KJ=  i  KnT  h /d)  t-K  - - - - 

JI  =  Jl  OH  ♦  J 

155 

=  <  f  THPT7P  (TtOH)  -  .1  >riCTCR7IL0W»l ,  JI  ,7317X1— TH7TT\j rirffRFFin  * - 

1  FACT  OP  (XLOH.JI  (KJ)  >/  ITNPTABIILOHJ-TMPTABFILON-fI)  I 

JI  -  JL  n  N  fj-j  “  - - - 

F«CIT(J)  *  (21) -MOLTHl  'FAC 

49  CuniiNi'E  '"  ""  "  . —  —  — ■ - - - - - - 

TFACT<r>  =  »F^CIK2)-FACITJ(l))/(ri.HD»(JLPH4i»-TLMD»(JLaMt7)> 

c 

50  C LN  1  I  Nl't  - - - - - - - - - 

CONFUTE  ATTENUATION  (OA/KN) 

PIIN/L  J  »  L  - - - - - - — - 

GMPAIN*  AlTM  RATFS.  ATTN.RATE.J?  * 

165 

l  ««K*  it-KP!  Imoapj  /  J T  T  r  flU  M  Z  lTT^Tir - - 

2 T FACT  <1  ‘ 0iKw^2)  K JJ  > 

HHJKN  - - 

FNO 

